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ABSTRACT 

The discovery of e Indi Ba, Bb, a binary brown dwarf system very close to the Sun, makes possible a concerted campaign to char- 
acterise the physical parameters of two T dwarfs. Recent observations suggest substellar atmospheric and evolutionary models may 
be inconsistent with observations, but there have been few conclusive tests to date. We therefore aim to characterise these benchmark 
brown dwarfs to place constraints on such models. We have obtained high angular resolution optical, near-infrared, and thermal- 
infrared imaging and medium-resolution (up to R~5 000) spectroscopy of s Indi Ba, Bb with the ESO VLT and present VRIzJHKL' M' 
broad-band photometry and 0.63-5.1 jim spectroscopy of the individual components. The photometry and spectroscopy of the two 
partially blended sources were extracted with a custom algorithm. Furthermore, we use deep AO-imaging to place upper limits on the 
(model-dependent) mass of any further system members. We derive luminosities of log L/Lq=-4.699 ± 0.017 and -5.232 ± 0.020 
for e Indi Ba, Bb, respectively, and using the dynamical system mass and COND03 evolutionary models predict a system age of 
3.7-4.3 Gyr, in excess of previous estimates and recent predictions from observations of these brown dwarfs. Moreover, the effective 
temperatures of 1352-1385 K and 976-101 1 K predicted from the COND03 evolutionary models, for b Indi Ba and Bb respectively, 
are in disagreement with those derived from the comparison of our data with the BT-Settl atmospheric models where we find ef- 
fective temperatures of 1300-1340 K and 880-940 K, for £ Indi Ba and Bb respectively, with surface gravities of log g=5.25 and 
5.50. Finally, we show that spectroscopically determined effective temperatures and surface gravities for ultra-cool dwarfs can lead to 
underestimated masses even where precise luminosity constraints are available.*** 

Key words. Stars: atmospheres - Stars: fundamental parameters - binaries: general - Stars: low-mass, brown dwarfs - Stars: late-type 
- Stars: individual: e Indi B 



1. Introduction 

The characterisation of low-mass stars and brown dwarfs is im- 
portant for studies of substellar and planetary atmospheres, the 
reliable application of low-mass evolutionary models, and the 
derivation of the full initial mass function. With over five hun- 
dred L and over one hundred T dwarfs now knowrQ statistical 
studies of global properties and detailed studies of the closest 
objects are now possible. 

Binary systems have an important role to play. They allow 
the determination of dynamical masses, provide a laboratory in 
which objects with the same age and chemical composition may 
be compared, and, where they have main-sequence companions, 
provide external constraints of metallicity and age which iso- 
lated objects lack, breaking the substellar mass-luminosity-age 
degeneracy. 



To fully constrain the evolutionary models of substellar ob- 
jects (e.g. Burrows et al.| 1997} Baraffe et al.||2003 Saumon & 



* Based on observations collected with the ESO VLT, Paranal, Chile 
under program 072.C-0689 
** Email: rob@astro.ex.ac.uk 
*** The full resolution spectra of both brown dwarfs are available in 
electronic form at the CDS via anonymous ftp to cdsarc.u-strasbg.fr 

1 http://dwarfarchives.org - the M, L, and T dwarf compendium main- 
tained by Chris Gelino, Davy Kirkpatrick, and Adam Burgasser. 



Marley 2008 ), it would be most useful to determine the bolo 
metric luminosity, radius, mass, and age of a range of such ob- 
jects. Bolometric luminosities can be determined from photo- 
metric and spectroscopic observations across a large wavelength 
range. Masses can be determined in systems where an orbit may 
be monitored, and finally, the age and metallicity can be inferred 
from better characterised stars in the same system. To constrain 
the atmospheric models of brown dwarfs, it is necessary to ac- 
quire high signal-to-noise spectra over as wide a wavelength 
range as possible, allowing robust estimates of the effective tem- 
perature and surface gravity to be made. 

The discovery of a distant companion (projected separation 
~1500AU) to the high proper-motion (~4.7 arcsec/yr) K4.5V 
star, e Indi, was reported by Scholz et al. (2003). One of 
our nearest neighbours, e Indi A has a well-constrained par- 
all ax from HIPPARCOS ( |Perryman & ESA|[L997l > as refined 
by |van Leeuwen (2007), putting the system at a distance of 



3.6224±0.0037pc. This was followed by the discovery of the 
companion's binary nature (McCaughrean et al. 2004). The 
proximity of s Indi Ba, Bb to the Earth means Ba is more than 
a magnitude brighter than any other known T dwarf, and allows 
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unprecedented, detailed spectroscopic studies of these important 
template objects. 

s Indi Ba, Bb are uniquely suited to provide key insights 
into the physics, chemistry, and evolution of substellar sources. 
Although there are a number of other T dwarfs in binary sys 
tems, such as the M4/T8.5 binary Wolf 940 (Burningham et al. 



2009) and the T5/T5.5 binary 2MASS 1534-2952 (Liu et al 



2008), s Indi Ba, Bb has a very well-determined distance, a 



main-sequence primary star with which to constrain age and 
metallicity, and a short enough orbit (nominally ~ 15 years, 
McCaughrean et al. 2004 1 such that the system and individ- 
ual dynamical masses can soon be determined (McCaughrean] 
eTaT1|2009| |Cardoso et aT1 |2009b| in prep.). They are also rel- 



atively bright, close enough, and sufficiently separated to allow 
detailed photometric and spectroscopic studies of both compo- 
nents. Importantly, these two objects roughly straddle the L to 
T transition (cf. Burgasser 2009 ) where the atmospheres of sub- 
stellar objects alter dramatically. The study of these two coeval 
objects on either side of the transition will help in understand- 
ing the processes effecting the change from cloudy to cloud-free 
atmospheres. Characterisation of this system allows the mass- 
luminosity-age relation at low masses and intermediate age to 
be tested, investigation of the atmospheric chemistry, including 
vertical up-mixing, and detailed investigation of the species in 
the atmosphere. 

To date, spectroscopic observations of T dwarfs have pre- 
dominantly been either at low-resolution (e.g., Burgasser et al. 
2002 Chiu et al.|2006 1, which allows spectral classification and 
overall spectral energy distribution modelling to determine lu- 
minosities, or high-resolution studies of relatively small wave- 
length regions to investigate gravity and effective temperature- 
sensitive features. For example, |McLean et al.| (|2003) presented 
near-IR spectra at a spectral resolution of R~2000 of objects 
spanning spectral types M6 to T8, and discussed broad changes 
in spectral morphology and dominant absorbers through the 
spectral sequence. This was complemented by R~20 000 /-band 
spectra presented in McLean et al. ( 2007 ) where many H2O and 
FeH features were identified and the progression of the /-band 
potassium doublet from M to T dwarfs charted. 

Previous studies of ultra-cool dwarfs attempting to constrain 
low-mass evolutionary models have been hampered by ambigu- 
ous ages, possible unresolved binarity, and the difficulty asso- 
ciated with acquiring observations of close, faint companions. 
For example, observations of AB Dor C have roused some con- 
troversy over the applicability of current low-mass evolutionary 
models, with the assumed age of the system being a major source 
of disagreement. |Close et al.| ( |2005[ ) determined a dynamical 
mass for AB DorC and, using an assumed age of 30-100 Myr, 
argued that evolutionary models predicted a higher luminosity 
than was observed. However, Luhma n~et al.| ( |2005| ) countered 
with an analysis based on an age of 75-150 Myr, finding no 
significant discrepancy between the observations and models. 
Nielsen et al.| (|2005|) further argued that using their slightly re- 



vised age of 70+30 Myr, the models still under-estimated the 



mass of this object. Again, this was disputed by Luhman & 
(|2006 1 after a re-reduction of the same data used by Close 



Potter 



et all 



2005| ), and then |Close et al.| ( |2QQ7] > concluded that, based 



on newly acquired spectra, there was no discrepancy between the 
observations and models. Despite this apparent rapprochement, 
the situation may nevertheless be further complicated by the sug- 
gestion of Marois et al. (2005} that AB Dor C may itself be an 
unresolved binary. More recently, Dupuy et al. (|2009) presented 
a dynamical mass for the binary L dwarf system HD 130948BC 
which, along with an age estimated from the rotation of the 



main-sequence parent star, suggests that the evolutionary models 
predict luminosities 2-3 times higher than those observed. 

LeggettetallpOOS] ) used 0.8-4.0 jum spectra at R~ 100^60 



and near- to mid-lR photometry of HN Peg B, a T2.5 dwarf com- 
panion to a nearby GOV star, to investigate physical properties 
including dust grain properties and vertical mixing. In the near- 
IR, the resolution was too low to study spectral lines in detail. 
However, by fitting the overall spectral morphology and making 
use of the longer-wavelength data, they were able to place impor- 
tant constraints on vertical mixing and sedimentation. Leggett 
|et al.| ( |2009[ ) also reported the physical properties of four T8- 
9 dwarfs from fitting observed near- to mid-IR spectral energy 
distributions with the atmospheric and evolutionary models of 
|Saumon & Marley| ( |2008[ ). They discussed the effects of vertical 
transport of CO and N2 and demonstrated the complementary 
effects of increasing metallicity and surface gravity. 

Reiners et al. ( |2007| ) analysed high-resolution (R~33 000) 



optical spectra of three L dwarfs and the combined s Indi Ba, 
Bb system, concluding that although some individual features 
are not well-matched by the model atmospheres and that sig- 
nificant differences remain for some molecular species and al- 



kali metal features, general features are reproduced. |Smith et al. 



( |2003[ ) also acquired high resolution (R~50 000) near-IR spectra 
of (only) e Indi Ba in the wavelength ranges 1 .553-1 .559 /mi and 
2.308-2.317 /im. These were fit with the unified cloud models of 
Tsqi] ( |2002) and effective temperatures of 1400 K and 1600 K 
were derived for the two spectral regions. Mid-IR spectroscopy 
of the unresolved s Indi Ba, Bb system was also acquired by 
Roellig et al. ( 2004[), who use evolution ary models along with 
the luminosities of McCaughrean et al. ( 2004 1 and an assumed 
age of 0.8-2.0 Gyr to derive effective temperatures and surface 
gravities and then compared composite spectral models to the 
observed spectrum. Their predictions were revised by Mainzer 
et al. ( 2007| ) who derived effective temperatures of 1210-1250K 
and 840 K, for s Indi Ba and Bb respectively, under the assump- 
tion of a system age of ~ 1 Gyr. 

Finally, |Kasper et H] ((2009 ) presented R-400 near-IR 
NACO/VLT spectroscopy of e Indi Ba, Bb which were com- 



pared to the evolutionary models of Burrows et al. ( 1997[ ) and 
the atmospheric models of Burrows et al. ( 2006| ). They derived 
effective temperatures of 1250-1300 K and 875-925 K, and sur- 
face gravities of log g (cms -1 ) 5.2-5.3 and 4.9-5.1, for e Indi Ba 
and Bb respectively, by comparing their observed spectra with 
their spectral models scaled using the distance and a radius pre- 
dicted by their evolutionary models. We will discuss these results 
further in Sect.[l0]in contrast to our new data. 

In this paper we present high signal-to-noise photometry 
from the V- to M'-band (0.5^.9 fim) and medium resolution 
spectroscopy from 0.6-5.1 //m of the individual components of 
the s Indi Ba, Bb system. In Sect. [2] we describe the observations 
and data reduction, including the routines employed to extract 
the partially-blended photometry and spectroscopy. We re-derive 
the spectral types of both objects according to the updated classi- 
fication scheme of Burgasser et al. (2006b ) in Sect.[3]and discuss 
constraints imposed by the parent main-sequence star in Sect. [4] 
We derive the luminosities of both sources in Sect.|5]and discuss 
the preliminary dynamical mass measurement of McCaughrean 
et al.| ( [2009 in prep.) in Sect.|6] Our observations are compared 
to evolutionary models in Sect.|7]and to atmospheric models in 
Sect. [8] We then put limits on the masses of lower-mass compan- 
ions in Sect. [9] and finally the predictions of evolutionary and at- 
mospheric models and previous determinations are compared in 
Sect.[T0] 
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2. Observations and Reduction 

s Indi Ba, Bb were observed with the ESO VLT using 



FORS2/UT1 (Appenzeller et al. 



1998) for optical photometry 



and spectroscopy, ISAAC/UT1 ( Moorwood et al. 1998[ ) for near- 
to thermal-IR photometry and spectroscopy , and NACO/UT4 
(Lenzen et al. 2003 1 Rousset et al. 2003| l for deep near-IR 
AO imaging. In all observations except those using NACO, the 
point-spread functions (PSFs) of the two sources were partially 
blended, even under excellent observing conditions with seeing 
always less than 0.7". 

2.1. optical photometry 

Broadband VRIz photometry was obtained on June 19 and July 
20 2004 (UT) using FORS2 (2 CCDs each 2048 x 4096 pix- 
els) in high resolution mode with 2x2 binning resulting in a 
plate-scale of0.125"pixeF 1 and a field-of-view of 4. 25' x4.25'. 
s Indi Ba and Bb were separated by ~0.84" under photomet- 
ric conditions with median seeing of 0.55" FWHM. Five im- 
ages dithered by 1' from a central position were obtained in 
each filter except the /?-band where twelve dithered images were 
taken. Individual exposure times were 500s, 60s, 20s, and 10s 
in the VRIz bands, respectively, giving total integration times 
of 42 min, 12 min, 100 s, and 50 s. Sky subtraction and flat- 
fielding were carried out with standard IRAF programs. As seen 
in Fig.[TJ both components of the binary are well-detected in the 
RIz bands, but s Indi Bb is only marginally detected in the V- 
band. We used the FORS2 Bessell V, Special R, Bessell /, and 
Gunn z broadband filters. Observations of the standard star fields 
PG 2213-006 and Mark- A ( |Landolt|1992] l were taken for photo- 
metric calibration which is discussed in detail in Sect.lCl 

The large field-of-view meant there were sufficient bright 
field stars with which to model the PSF and so DAOPHOT/IRAF 
PSF-fitting was employed to extract the individual fluxes of the 
two brown dwarfs. For the RIz bands, each of the images were 
fit separately to allow a determination of the accuracy of the pro- 
file fitting which is included in the uncertainties of the derived 
magnitudes. In the V-band, s Indi Bb was only marginally de- 
tected, so we were unable to fit PSFs to the two components 
of the binary. To extract the photometry of both sources, we 
used a small aperture to measure the flux of the brighter source 
ensuring there was no appreciable contaminant flux from the 
fainter source. We then extracted the photometry of the com- 
bined source with a larger, circular aperture and, using the curve 
of growth of brighter stars in the field, derived the excess flux 
due to s Indi Bb and thus an upper limit on the V-band flux. 

The measured flux ratio and the central wavelength and 
width for each of the observed filters is listed in Table [T] while 
Table[2]lists the derived photometry. 

2.2. near- and thermal-IR photometry 

The ISAAC imager was used on November 5 and 1 1 2003 (UT) 
to obtain photometry of s Indi Ba, Bb in the ISAAC JHK S LM NB 
filters. ISAAC was used with the ALADDIN array (1024 x 1024 
pixels) in long-wavelength imaging modes LWI3 and LWI4 for 
the JHKs and LMnb imaging, respectively, resulting in plate- 
scales of 0.148" pixel -1 and 0.071" pixel -1 and fields-of-view 
of 151x151" and 73x73". 

s Indi Ba and Bb were separated by ~0.77" at this epoch and 
were observed under photometric conditions with typical seeing 
of ~0.45" FWHM. The JHKs images were taken at three posi- 
tions dithered by ~14" with total integration times of 3.5 min in 
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Fig. 1. From left to right and top to bottom, FORS2 VRIz and 
ISAAC JHKsLMnb images of s Indi Ba, Bb. Each image is a 
4" x 4" sub-section of the full image. North is up, East left. The 
object to the south-east of the V-band image is a faint galaxy not 
seen in the other pass-bands. The L-band profiles are seen to be 
slightly elliptical. The V and M NB -band images are the stacked 
images of all observations in those pass-bands. 

Table 1. The flux ratio (Ba/Bb) of e Indi Ba, Bb in the observed 
FORS2 and ISAAC filters and the effective wavelength and half- 
power width of those filters. See Sect.|C]for a description of pho- 
tometric calibration and transformation to standard systems. All 
except the V-band ratio were measured using PSF-fitting (see 
Sect.|2~T]i. 



Observed Passband 


Ratio (Ba/Bb) 


/irff (jjm) 


A/i (jim) 


V 


9.82 ± 


0.66 


0.554 


0.112 


R 


4.81 + 


0.07 


0.655 


0.165 


I 


5.04 ± 


0.05 


0.768 


0.138 


z 


3.85 ± 


0.03 


0.910 


0.131 


.1 


2.28 ± 


0.02 


1.25 


0.29 


H 


5.27 ± 


0.06 


1.65 


0.30 


K s 


7.12 + 


0.02 


2.16 


0.27 


L 


4.50 ± 


0.16 


3.78 


0.58 


M m 


1.35 ± 


0.15 


4.66 


0.10 



each of the three filters. The three offset images were combined 
to remove the sky flux and then flat-fielded with standard IRAF 
programs. The L and M NB images were taken in chop-nod mode 
with a throw of ~20" and total exposure times of 2 min in the L- 
band, and 4 min in the MNB-band. The images were flat-fielded 
and the half-cycle frames subtracted in the standard manner pro- 
ducing three sky-subtracted images at different positions on the 
array. 

The ISAAC field-of-view was found to be mostly empty, and 
most importantly with no stars bright enough to fit a model PSF. 
Therefore, a custom profile-fitting routine was implemented to 
determine the flux ratio of the partially-blended objects (see 
Sect. [A}. The total flux of the two objects was then measured 
by aperture photometry. The centre and size of a circular aper- 
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Table 2. The derived apparent magnitudes for e Indi Ba, Bb. 
The optical magnitudes are in the FORS2 system as described 
in Sect.[Cj and the near-IR magnitudes have been transformed 
into the MKO JHKL'M' filters via synthetic photometry. The 
uncertainties for both objects are dominated by the uncertain- 
ties on the measured flux ratio and the standard star magnitudes. 
Absolute magnitudes may be derived using the distance modulus 
of -2.205. 



Passband 


Ba 




Bb 


V 


24.12 


± 0.03 


> 26.60± 0.05 


R 


20.65 


±0.01 


22.35 ± 0.02 


I 


17.15 


±0.02 


18.91 ± 0.02 


z 


15.07 


±0.02 


16.53 ± 0.02 


.1 


12.20 


±0.03 


12.96 ± 0.03 


H 


11.60 


±0.02 


13.40 ±0.03 


K 


11.42 


±0.02 


13.64 ±0.02 


U 


9.71 


±0.06 


11.33 ±0.06 


M' 


10.67 


±0.23 


11.04 ±0.23 



ture were chosen so as to minimise the sky noise contribution 
while ensuring the total flux was not biased toward either object 
by including substantially more of the profile wings of one ob- 
ject than the other. The photometric calibration is discussed in 
detail in Sect. [C] and the derived near-IR photometry is listed in 
Table |2] 

2.3. photometric variability 



Koen| ( |2005| ) detected 7 c -band variability in optical photometry 



of the combined e Indi Ba, Bb system on time-scales of hours, 
reporting a linear rise in the Zc-band of ~0. 16 m over the course 
of 3.6 hours. Table[3]shows the /-band magnitudes for the com- 
bined system from different studies. For comparison, our obser- 
vations have been combined to give the magnitude of the unre- 
solved system. However, a direct comparison is hindered by the 
different response functions employed in these measurements. 
We find that a shift in filter response of 100 A is sufficient to ex- 
plain the spread in the /-band magnitudes due to the steep rise in 
T dwarf optical spectra, and so may shield any intrinsic variabil- 
ity. 

Evidence for variability of s Indi Ba, Bb of ~0.05 m was 
also found by Koe n et aL] ( |2005| l in the near-IR. Table |4] shows 
the 2MASS JHK$ magnitudes extracted from our flux cal- 
ibrated spectra (see Sect.|C]) and the 2MASS magnitudes of 
|McCaughrean et al.| ( |2004| i and |Kasper et aL] ( |20T)9l >. These 
mostly agree within the stated uncertainties with the exception 
of the //-band magnitude of e Indi Bb. However, the photome- 
try of McCaughrea n et al.| and [Kasper et akj was acquired in the 
VLT/NACO system with an //-band filter which extends into the 
region of high telluric water absorption, possibly accounting for 
the differing results. The uncertainties on our near-IR photome- 
try are of similar magnitude to the proposed variability. 

We plan to use the near-IR photometry (JHKs) obtained as 
part of the astrometric monitoring of this binary to further in- 
vestigate any variability over the monitoring epochs (-monthly 
from August 2003 to present) to probe longer-term variability 
than that detected by Koen et al. ( 2005| >. 



2.4. AO deep companion search 

NAOS/CONICA (NACO) was used on November 7 2003 (UT) 
to obtain deep adaptive-optics (AO) //-band imaging of the field 



Table 3. The /-band magnitudes adopted for the combined e Indi 
Ba, Bb system by different studies. We have used the FORS2 
Bessell /-band filter which is similar to the DENIS /-band, while 
Koen (2005) used the Cousins /-band. 



Date 



e Indi Ba, Bb Reference 



1997.771 
1999.666 
2000.781 
2004.462 
2004.546 



16.59±0.10 
16.77±0.10 
16.90±0.12 
16.70±0.04 
16.95±0.03 



SSS-UK 
SSS-UK 
DENIS 



Koen 
TRTs 



l2005|> 
work 



Table 4. The 2MASS magnitudes adopted for s Indi Ba , Bb in 
different studies. McC04 refers to |McCaughrean et al.| ( |2004} 
for which we have assumed the uncertainties on the photometry 
to be equal to the uncertainty on the 2MASS photometry of the 



unresolved system. K09 is the study of Kasper et al. (2009 



Passband 


Ba 












McC04 




K09 




This work 


2MASS J 
2MASS H 
2MASS K s 


12.29 ± 
11.51 ± 
11.35 ± 


0.02 
0.02 
0.02 


12.33 ± 
11.54 ± 
11.37 ± 


0.02 
0.02 
0.02 


12.29 ± 0.03 
11.50 ±0.03 
11.38 ±0.03 




Passband 


Bb 












McC04 




K09 




This work 


2MASS J 
2MASS H 
2MASS K s 


13.23 ± 
13.27 ± 
13.53 ± 


0.02 
0.02 
0.02 


13.19 ± 
13.36 ± 
13.51 ± 


0.03 
0.03 
0.02 


13.23 ± 0.03 
13.19 ±0.03 
13.49 ± 0.03 



around s Indi Ba, Bb. The N90C10 dichroic was used to send 
10% of the source flux to the science camera and 90% to the 
IR wavefront sensor (IR WFS). The median natural seeing con- 
ditions were 0.58" FWHM at the time of our observations and 
with the AO correction we obtained a FWHM of 0.12" in the 
final combined image. 

The S27 camera was used with a plate-scale of 
0.027" pixeF 1 resulting in a field-of-view of 27.7" x27.7". We 
obtained a total of 72 individual images with exposure times of 
135 s, giving a total integration time of 162 min. The 9-point 
dither pattern used resulted in a continuous radial coverage of 
10.94" measured from the position of s Indi Ba, corresponding 
to 39.6 AU at 3.622 pc. Sky subtraction, flat-fielding, image 
alignment and stacking were standard. 

As the two sources were well-separated, we were able to em- 
ploy DAOPHOT/IRAF to fit a PSF to s Indi Ba and use this 
as a model to subtract Bb, iteratively fitting the PSF and sub- 
tracting until Bb was well-removed. This image of Ba alone was 
then used to add scaled and offset objects into the original im- 
age to investigate our detection limits. We find no sources from 
7 AU (71 pixels, 1.93") out to the image edge, corresponding to 
39.7 AU, with flux greater than 0.1% of e Indi Ba, correspond- 
ing to a peak pixel flux 5cr above the background, or a source 
//-band magnitude of 19. l m which would have been detectable 
in our deep image. 

Closer to s Indi Ba, Bb, the noise increases due to the 
Poisson noise from the source flux and so the brightness lim- 
its on any companions are higher. The pixelation of the profile 
also acts to suppress the visibility of close companions. We find 
no sources with a flux greater than 2% of the s Indi Ba flux 
(Z/=15.8 m ) down to 0.83 AU (8.5 pixels, 0.23") from Ba or Bb, 
and down to 0.4 AU (4.3 pixels, 0. 12") from each object, we can 
discount any sources with more than 10% of the s Indi Ba flux 
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(//=14.1 m ). In Sect.[9]we use these flux limits to derive the mass 
limits on any possible companion. 

2.5. optical spectroscopy 

FORS2 (2 CCDs each with 2048 x 4096 pixels) was used on 
June 16 2004 (UT) to obtain optical (0.63-1. 07 pm) spec- 
troscopy of s Indi Ba, Bb in long-slit mode with a 30" dither 
along the slit. We used a 0.5" wide slit, the HR collimator, and 
2x2 binning mode, resulting in a plate-scale of 0.125" pixel -1 . 
The 600 RI and 600 Z grisms were used to obtain the full 0.63- 
1.07 pm spectrum, yielding resolutions of R~ 1000 at 6780 A and 
R-2000 at 9700 A, while the median seeing of 0.47" FWHM al- 
lowed the spectra of the two objects to be resolved at this epoch 
when the separation was ~0.84". Total integration times of 80 
min with the 600 RI grism and 38 min with the 600 Z grism 
were obtained by co-adding 6x800 s and 5x460 s individual ex- 
posures, respectively. 

The partially-blended spectra were extracted using a profile- 
fitting routine as described in Sect.|B] Wavelength calibration 
used arc spectra for the 600 RI grism and skylines for the 600 Z 
grism, with excellent agreement in the crossover region. Our 
spectra have vacuum wavelengths. There was no visible dif- 
ference between the two resolutions in the crossover region 
(0.75-0.85 pm), so they were combined without smoothing. 
Flat-fielding was carried out with lamp flats. Observations were 
also made of the white dwarf LTT9491 for relative flux calibra- 
tion. However, these data were taken with a wider slit to the s 
Indi Ba, Bb spectra and so could not be used for removal of tel- 
luric absorption. A spline was fit to the observed white dwarf 
spectrum which was used along with the known spectrum of 
LTT9491 ( |Okel[T990l [Ramuy et al.||1994} to model the wave- 
length dependence of the flat-field and detector response and so 
remove this from our observed target spectra. 

We scaled the NSO Kitt Peak atmospheric transmission 
~d2003l 



spectrum of |Hinkle et al. ( |2003[ ) to construct a grid of atmo- 
spheric spectra with varying absorption strengths (assuming all 
absorbing species vary similarly with airmass) to model the at- 
mosphere above Paranal at the time of our observations. These 
were Gaussian smoothed to the resolution of our observed spec- 
tra and the atmospheric model scaled to maximise removal of 
known features. The telluric features appear to be fully removed 
in the sections of the spectrum where we expect few features 
intrinsic to T dwarfs, which leads us to believe that we have ac- 
ceptable telluric removal even across the 9300-9800 A region 
where there are blended telluric and intrinsic features. 

The final optical spectra have a peak signal-to-noise per pixel 
of ~300 and ~200 for e Indi Ba and Bb, respectively, falling to 
-35 and -20 at 7000A. The 0.6-5.1 /jm spectra of e Indi Ba, 
and Bb are shown in Figs.[3]and|4]in both linear and logarithmic 
units, and the full resolution spectra are shown in Figs. [3] [6] and 
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2.6. near-IR spectroscopy 

ISAAC and its HAWAII array (1024 x 1024 pixels) with a plate- 
scale of 0.146"pixeF' was used on November 8 2003 (UT) in 
short-wavelength medium-resolution (SWS1-MR) mode to ob- 
tain near-IR spectroscopy from 0.9-2.5 pm of s Indi Ba, Bb with 
a resolution of R~5000. By turning the instrument rotator, both 
sources (at a separation of ~0.77" at this epoch) were simultane- 
ously placed on the 0.6" wide slit with median seeing of 0.50" 
FWHM. Three 60s exposures dithered by 20" along the 120" 




1.5 2 
Wavelength (|im) 



2.5 




1.5 2 
Wavelength (fxm) 

Fig. 2. Our spectra of s Indi Ba and Bb (red lines), smoothed to 
a resolution of 17A FWHM and median filtered in the regions of 
high telluric absorption, compared to the lower resolution spec- 
troscopy of |Kasper~ ;t al. (2009) (green lines). It is seen that the 
overall match of the absolute scale is reasonable. Here we have 
absolutely flux-calibrated using our ISAAC photometric obser- 
vations, while |Kasper et alT| ( |2009[ ) use the 2MASS point source 
catalogue magnitudes. 



long slit were taken in each of twenty-one wavelength regions 
to cover the entire 0.9-2.5 pm range. We observed three spectral 
regions in the 0.98-1. 10pm domain each 0.046 pm wide, and 
six spectral regions in each of the 1.10-1. 40 pm, 1.40-1. 82 pm, 
and 1.82-2. 50 pm spectral domains each covering a range of 
0.059 pm, 0.079 pm, and 0.122 pm, respectively, with a mini- 
mum cross-over between regions of 0.004 pm allowing full 0.9- 
2.5 pm spectra to be compiled with no gaps in coverage. 

Observations were also made of the nearby G2 dwarf 
HD209552 in each of the wavelength regions interspersed with 
the target observations in order to flux calibrate and correct for 
telluric absorption as first described by |Maiolino et al. ( |1996 1. 
The airmass difference between the observations of s Indi Ba, Bb 
and the telluric standard were in the range 0.07-0.17. Tungsten- 
illuminated spectral flats were taken in the same configuration 
at the end of the night and wavelength calibration was achieved 
through use of OH skylines for all but four of the wavelength re- 
gions which had too few (known) skylines for a reasonable fit. In 
these cases, XeAr arc spectra were used and cross-over regions 
ensured a good match throughout. 
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Fig. 3. Full 0.6-5.1 //m spectrum of s Indi Ba and Bb (upper and lower lines respectively). Here the full resolution optical and 
near-IR spectra have been smoothed to 17 A FWHM to allow inspection of the broad features without being dominated by the many 
fine features seen in the full resolution spectra (see Figs. BUT}. 
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Fig. 4. The same as Fig.[i]but with the flux on a logarithmic scale. 



For each of the wavelength regions, the three dithered spec- 
tral images were combined to subtract the sky background. This 
process leaves residuals of the order of a few percent of the ob- 
served counts in the spectrum caused by temporal variations 
in the strength of the sky emission, negligible for all but the 
strongest sky lines. The residual sky was further subtracted by 
taking apertures above and below the spectra and subtracting 
the sigma-clipped mean sky level. The images were then di- 
vided by the spectral dome flats. The partially-blended spectra 
were extracted in the same manner as the optical (see Sect.[B|. 
The extracted spectra of s Indi Ba and Bb were divided by 
the standard star spectrum and corrected using the solar spec- 
trum of Walla ce~et al.| ( |1996| l smoothed to the resolution of 
our observations to remove the wavelength dependence of the 
detector and the flat-field. In the regions of very high telluric 
absorption between the / and H and H and K bands where 
there were gaps in the observed solar spectrum (13517-14032 A 
and 180 24-193 17 A), w e substituted the Kurucz IRRADIANCE 
modej^] (Kurucz 2005 ) to avoid any breaks in the coverage of 



our spectrum. These regions are generally omitted from pub- 



lished spectra, but since we have high-enough signal-to-noise 
and spectral resolution, we retain them to allow a comparison of 
the continuum flux in these regions to model predictions. 

The known spectral type of the standard star, response of the 
filter systems used, and magnitudes of the standard stars allowed 
us to flux calibrate a template spectrum for each standard star in 
each of the JHKL'M' passbands. We then employed the ratio 
of the observed fluxes of the targets and standard stars in the 
ISAAC system, and the known filter responses to flux calibrate 
our observed e Indi Ba and Bb spectra. This was done separately 
for each passband as each region is bounded by high telluric ab- 
sorption, but the optical spectra were flux calibrated by scaling 
to match the /-band spectrum in the 9782-9999 A cross-over 
region. Figure|2] shows a comparison of our final near-IR spec 
tra for both brown dwarfs with those from Kasper et al. ([2009 



http://kurucz.harvard.edu/sun/irradiance/solarirr.tab 



No scaling has been applied to these spectra which, nonetheless, 
have similar absolute flux calibrations. While there are some 
small discrepancies between the spectra, the overall match sug- 
gests that no large systematic error has been introduced by the 
combination of many short wavelength regions to produce our 
full near-IR spectrum. The final spectra have a signal-to-noise 
of -80/60 (Ba/Bb) per pixel at the /-band peak, 100/70 at the 
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//-band peak, and ~ 100/50 at the peak of the /f-band spectrum. 
The full resolution observed spectra of both objects are available 
at CDS and on the author's web-pagesj^] 

2. 7. thermal-IR spectroscopy 

ISAAC and its ALADDIN array (1024 x 1024 pixels) with a 
plate-scale of 0.146" pixel -1 was used on November 6-7 2003 
(UT) in long-wavelength, low-resolution (LWS3-LR) mode to 
obtain R~600 and R~500 spectroscopy of s Indi Ba, Bb in the 
spectral ranges 2.86^4.19//m (L-band) and 4.53-5.08 yum (M- 
band), respectively, using slit widths of 0.6" and 1.0". The M- 
band spectra had to be binned before they could be extracted giv- 
ing a final resolution of R~220 at 4.75 fim. Both sources were 
placed on the 120" long slit and chop-nod mode was used re- 
sulting in a total on-source time of 30 min and 35 min for the L- 
and M-bands, respectively, with median seeing of 0.54" FWHM. 
The half-cycle frames were subtracted in the standard manner 
producing three sky-subtracted spectra at different positions on 
the array. Observations were also made of the solar-type star 
HD2 10272 in the same manner to allow flux calibration and cor- 
rection for telluric absorption as for the 0.9-2.5 fim spectra. The 
allace et al. ( 1996| > solar spectrum was incomplete with a gap 
in the range 4. 17-4.55 /im which affected the last ~0.02yum of 
the L-band spectrum and the first ~0.02 [im of the M-band spec- 
trum. Again, the gap was filled by the Kurucz IRRADIANCE 
model. Tungsten flats were taken with the same slits at the end 
of the night and wavelength calibration was achieved through 
use of XeAr arc spectra. Section|B] explains in detail how the 
partially-blended spectra were extracted. The final spectra have 
a peak signal-to-noise of ~40/15 (Ba/Bb) per pixel in the L-band 
and ~10 per pixel for both objects in the M-band. 

3. Spectral Classification 

3.1. near-IR spectral classification 



McCaughrean et al.| (|20Q4| > used the |Geballe et ah] ( |2002| i and 
Burgasser et al. (2002) classification indices to provide spectral 
classifications for e Indi Ba, Bb based on their //-band spec- 
tra. They arrived at spectral types of Tl and T6 for Ba and Bb, 
respectively, by employing both of the Burgasser et al. (2002) 
//-band indices and the CH4 index of |Geballe et al.| ( |2002| ) The 
H2O index of Geb alle et al.| indicated TO and T4 and was ex- 
cluded based on previous spurious spectral classification of Gl 
229 B. 

Given the greatly improved data here, we employ both meth- 



ods of near-IR spectral index classification of Burgasser et al. 
(2006b), namely direct comparison of spectral indices with stan- 
dard indices and against index ranges defined for each subtype. 
In addition, we plot our s Indi Ba, Bb spectra alongside the spec- 
tra of the defined standards to allow direct morphological com- 
parison. 

Since the precise value of spectral indices is dependent on 
the spectral resolution, we calculate the indices with our spectra 
smoothed to resolutions of R=150 and R=500 equivalent to the 
resolutions for the standard stars as observed with IRTF/SPEX 
(at R-150) and UKIRT/CGS4 (at R-500). Table [5] shows the 
five spectral indices measured from the different resolution spec- 
tra for both s Indi Ba and Bb. We also list the values for the 
unsmoothed spectra to quantify the effect of comparing indices 
for spectra at different resolutions. Spectral types inferred from 
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Fig. 8. Smoothed spectrum (30 A FWHM) of s Indi Ba (thick 
lines) and spect ra of t he TO (SDSS 120747.17+024424.8; 
|Looper et al.| |2007[), Tl (alternate standard SDSS 



015141.69+124429.6; jBurgasser et al. 2004), and T2 (SDSS 
J125453.90-012247.4; |Burgasser et al.|2004| ) spectral standards. 
The three panels show the JHK bands omitting the lower 
signal-to-noise H2O absorption regions between the bands. 
Spectra are normalised to unity at 1.27/im in the /-band, at 
1.58 yum in the //-band, and at 2.10/mi in the Zf-band and offset 
to show the e Indi Ba spectrum against each of the spectral 
standards. 
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Fig. 9. Smoothed spectrum (30 A FWHM) of e Indi Bb (thick 
lines) and spectra of the T5 (2MASS 15031961+2525196 



(Burgasser et al. |2004|), T6 (SDSS 162414.37+002915.6; 


Burgasser et al. 


2006a), and T7 (2MASS 0727182+171001; 


Burgasser et al. 


2006a) spectral standards. Plotted regions and 



normalisation are as in Fig. [8] 



each index are shown in parentheses. The range of index values 



http://www.astro.ex.ac.uk/people/rob/Research 



for each subtype of the T class (method 2 of Burga sser et al 
|2006b) has only been defined for the R-500 CGS4 spectra of 
the spectral standard stars. 

Figure HJshows the spectrum of s Indi Ba smoothed to R= 150 
to match the resolution of the spectra of the standards and nor- 
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1.24 1.26 
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Fig. 5. The full resolution spectra of s Indi Ba (top line) and Bb for 0.63-1.3 fim. The FORS2 optical data is joined to the higher- 
resolution ISAAC near-IR data at 0.978 fim where we start to see even finer spectral features due to the higher spectral resolution 
and high signal-to-noise. Comparison of the two spectra reveals much of the apparent "noise" to be real spectral features present in 
both objects. The resolution is 6.8 A FWHM up to A = 0.79 fim, 6.5 A FWHM for A =0.79-0.978 fim, 1.8 A FWHM for A =0.978- 
1.1 fim, and 2.4 A FWHM for A =1.1-1.3 fim. Key spectral features are labelled. See McLean et al. (2003 i for a detailed discussion 
of features identified in brown dwarf spectra. 
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Fig. 6. Same as Fig. \^ but resolution is 2.4 A FWHM for A =1.3-1 A /mx, 3.2 A FWHM for A = 1.4-1. 82 //m, and 4.9 A FWHM 
for A =1.82-2.0//m. The regions 1.40-1.41 fim and 1.82-1. 92 fim have many artefacts due to high levels of telluric contamination. 
However, the high signal-to-noise and spectral resolution of these data allow us to extract the object spectra between the water lines, 
and so even in these regions we have a measure of the continuum flux. 
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Wavelength (|J.m) 

Fig. 7. Same as Fig. [5] but the observed spectra have gaps 2.50-2.86 jum and 4. 19^4.53 yum between the near-IR and the L- and 
M-bands. The resolution is 4.9 A FWHM for A =2.0-2.5 jtfln, 60 A FWHM for A =2.86-4. 19yum, and 216 A FWHM for A =4.5- 
5.1 fim. For clarity, the flux of e Indi Bb has been increased by a factor of 2 in the range 2.0-2.2 fim and by a factor 5 in the range 
2.2-2.5 fim. The telluric features in the M'-band (marked with ©) are due to very bright sky emission which dominates the source 
signal. The artefacts seen in the region 2.43-2.50 yt/m are again due to high levels of telluric contamination. 



Robert R. King et al.: s Indi Ba, Bb: a detailed study 



11 



Table 5. Spectral indices and classifications for s Indi Ba, Bb. The spectra have been smoothed to two resolutions to emulate the 
observations of the spectral standard stars defined by Burgasser et al. (2006b). The indices derived from our unsmoothed data have 



been included to quantify the uncertainty when measuring indices from spectra at different resolutions. Ml refers to method 1 of 
Burgasser et al. where the measured spectral indices are compared to that of each of the standard stars, and M2 refers to their method 
2 where the measured spectral indices are compared to a range of values for each spectral type. Direct classification compares the 
overall spectral shape to the standard star spectra. The spectral type derived from each index is shown in parentheses after each 
measurement with only Ml being defined for a resolution of R=500. The average spectral type from each method is shown on the 
right-hand side. Only direct classification is possible for our unsmoothed spectra. The range of the H2O-/ index is ill-defined for 
the earliest T dwarfs and is therefore omitted for s Indi Ba. 



Source 


H 2 0-J 


CH 4 -y 


H 2 0-H CH A -H 


CH 4 -K 


Direct 


Ml 


M2 


Spectra smoothed to R = 150 


e Indi Ba 


0.650 (T0.0) 


0.656 (T0.0) 


0.592 (T0.0) 0.873 (T2.0) 


0.587 (T2.0) 


T1.0 


T1.0 




e Indi Bb 


0.172 (T6.0) 


0.305 (T6.5) 


0.308 (T5.5) 0.295 (T6.0) 


0.178 (T5.5) 


T6.0 


T6.0 




Spectra smoothed to R = 500 


e Indi Ba 


0.652 (.../...) 


0.657 (T1.0/T1.5) 


0.589 (T1.0/T0.5) 0.872 (T2.0/T2.0) 


0.588 (T2.0/T2.0) 


T1.0 


T1.5 


T1.5 


£ Indi Bb 


0.171 (T6.0/T6.0) 


0.302 (T6.0/T6.0) 


0.303 (T6.0/T6.0) 0.290 (T6.0/T6.0) 


0.180 (T5.5/T5.5) 


T6.0 


T6.0 


T6.0 


Unsmoothed spectra 


£ Indi Ba 


0.645 


0.660 


0.589 0.870 


0.589 


T1.0 






£ Indi Bb 


0.171 


0.307 


0.299 0.287 


0.178 


T6.0 







malised within each of the JHK spectral regions, omitting the 
low signal-to-noise H2O absorption regions between the bands. 
The spectrum is plotted along with the spectra of the TO, Tl, 



and T2 dwarf standards observed with IRTF/SPEX ( |Looper et al. 
|2007| [Burgasser et al.|200 4). The Tl standard provides the best 
overall match, however there are some obvious discrepancies. 
The /-band spectra are normalised to one at 1.27/im and al- 
though the Tl standard fits this peak and the 1.15 /vm CH4/H2O 
absorption band well, the relative strength of the 1 . 1 //m peak is 
not well-matched. The //-band is again best matched by the Tl 
standard with however some deficit in the flux of e Indi Ba in 
the 1.62-1. 74 /mi region which seems to indicate slightly differ- 
ing levels of CH4 absorption. The Zf-band is better fit by the T2 
standard, however all three standard spectra struggle to match 
the region beyond 2.3 //m. 

Similarly, Figure [9] shows the smoothed spectrum of e Indi 
Bb along with the spectra of the T5, T6, and T7 dwarf standards 
observed with IRTF/SPEX ( |Burgasser eiliLpXHl [2006a ). The 
T6 standard provides the best overall match, however there are 
again some obvious discrepancies. The /-band spectra are nor- 
malised at 1.27 /mi and although the T6 standard fits this peak 
and the CH4/H2O absorption band well, the relative strength of 
the 1.1 yum peak is better matched by the T5 standard (although 
the absorption dip then no longer matches). The //-band is very 
well-matched by the T6 standard, with the T5 and T7 standards, 
respectively, over- and under-estimating the flux in the 1.62- 
1 .74 //m region. The Zf-band is also best fit by the T6 standard, 
although again all three standard spectra differ slightly beyond 
2.3 /mi. 

The discrepancies between the relative strength of the 1 . 1 /mi 
and 1 .25 /mi peaks between s Indi Ba, Bb and the spectral stan- 
dards shown in the left-most panels of Figs. [8] and [9] could possi- 
bly be explained by the effects of different metallicities, surface 
gravities, or cloud cover. This may further explain the appar- 
ent excess flux in e Indi Ba beyond 2.1 /mi relative to the Tl 
standard. We discuss the effects of surface gravity and slightly 
sub-solar metallicity of model spectra in Sect. |8.2| 

Finally, the derived near-IR spectral types for both objects 
from each of the three methods are summarised in Table [5] We 
adopt final near-IR spectral types of T1-T1.5 and T6 for e Indi 
Ba and Bb, respectively. The uncertainty on the spectral type of 
s Indi Ba is not due to poor signal-to-noise in our data nor the 
spectra of the standards, but rather is presumably due to second 



Table 6. Optical spectral indices and classifications for s Indi 
Ba, Bb. The spectra have not been smoothed to derive these in- 
dices. Direct classification compares the overall spectral shape 
with the standard star spectra. The spectral type derived from 
each index is shown in parentheses after each measurement. The 
spectral range of the FeH(B) index is limited to types later than 
T2 and the range of the colour-e index to types earlier than T2. 
The approximate classification from direct comparison for Ba is 
due to the lack of a TO or Tl optical spectral standard. 



Index 


£ Indi Ba 


£ Indi Bb 


CsI(A) 


2.147 (T2.0±1.0) 


1.738 (T6.0±1.0) 


CrH(A) / H 2 


0.906 (T0.0±1.0) 


0.452 (T6.0±1.0) 


FeH(B) 


1.237 (<T4?) 


1.109 (T6.5±0.5) 


Colour-e 


3.323 (T0.0±1.0) 


4.135 (> T2) 


Mean 


T0.5 


T6.0 


Direct 


<T2.0 


T6.0 



order spectral variations due to differences in metallicity and sur- 
face gravity. 



3.2. optical spectral classification 

We have also derived spectral types from our optical spectra 
based on the optical T dwarf classification scheme of Burgasser 
|et al.| ( [2003] l. The spectral indices are all defined for wavelengths 
above 0.90 //m and so only use the higher signal-to-noise portion 
of the standard star spectra. Table [6] shows the value of the four 
indices calculated from our full resolution spectra, along with 
the derived mean spectral type and the classification from direct 



comparison to the spectral standards defined by Burgasser et al. 
We do not smooth the spectra for this comparison. The direct 
spectral comparison is shown in Figs. [TO] and [TT| 

The optical classification broadly agrees with the near-IR 
classification for these two T dwarfs with spectral types of T0- 
T2 for s Indi Ba and T6.0-T6.5 for s Indi Bb. It is clear that our 
high signal-to-noise optical spectra will be useful as standards 
for future comparisons of optical T dwarf spectra. 
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Fig. 1 0. Optical spectrum of s Indi Ba and the T2 (SDSS 1254- 
0122 (Burgasser et al.|2003 i) optical spectral standard. There is 
no object defined as either the TO or Tl spectral standard in the 
optical. The s Indi Ba spectrum is normalised to unity at 9000 A 
and the T2 standard spectrum is offset by 1 dex for clarity. 




6500 



7000 



7500 



8000 8500 
Wavelength (A) 



9000 



9500 



10000 



Fig. 1 1. Optical spectrum of e Indi Bb and the T5 (2M ASS 0559 
1404;|Burgasser et al.|2003) , T6 (SDSS 1624+0029; |Burgasser 
et al.|2000| >, and T8 (2MASS 0415-0935; [Burgasser et al.|2003f 
optical spectral standards. There is no object defined as the T7 
spectral standard in the optical. Spectra are normalised at 9000 A 
and offset in steps of 1 dex for clarity. 



4. Constraints from e Indi A 

The binary s Indi Ba, Bb is a common proper motion compan- 
ion to the K4.5V primary, s Indi A, which allows us to con- 
strain some of the fundamental properties of the s Indi system. 
Most obviously, since the parallax of the primary was measured 
by HIPPARCOS ( |Perryman & ESA|[T997] l, we have an accu- 



rately known distance. This was updated in the reanalysis of 
|Leeuwen] ( |2007] > to yield 3.6224±0.0037pc. The uncertainty on 
the distance to s Indi Ba, Bb is somewhat larger: although we 
know the projected separation of the primary and brown dwarf 
binary to be 0.007 pc, we do not know the line-of-sight separa- 
tion. This added uncertainty should be resolved by the ongoing 
absolute astrometric monitoring which along with the mass ratio 
of the s Indi Ba, Bb system, allows determination of the parallax 
of the brown dwarf binary. 



Similarly, determinations of the metallicity of the primary 
can be applied to the brown dwarf companions assuming them 
to be co-eval and born from the same molecular material. The 
metallicity of s Indi A has been studied by a number of authors 
who arrive at somewhat differing results due to the chosen spec- 
tral lines and the different models with which the spectra are fit 
to derive the abundances. |Abia et al.| ([T988) derived a metallic- 
ity of [Fe/H]=-0.23, while the Geneva group reported metallic- 
ities in the range [Fe/H]=-0.2-+0.06 (Sant os et al.|2001| [2004) 1 
with their most recent determination of [Fe/H]=-0.2 reported 
by Sousa et al. (2008 1. These results suggest that s Indi A, and 



by association e Indi Ba, Bb, appear to have slightly sub-solar 
metallicity. In Sect. 8.2 we show atmospheric models with both 
[M/H]=0.0 and -0.2 fit to our spectroscopic observations. Note 
however, that [M/H] refers to the global metallicity and also de- 
pends on the abundance of ar-elements (cf. Ferraro et al. 1999 ), 
so is not necessarily equivalent to [Fe/H]. 

In previous analyses, the age estimate of 0.8-2.0 Gyr for 
s Indi A from |Lachaume et al.| ( |1999| l was used to make 
predictions from evolutionary models. Using this age range 
|McCa ughrean et al. (2004) derived model masses of 47+10 and 
28+7 Mj up . As we will discuss, we now have direct determina- 
tions of the luminosity of both brown dwarfs and of the system 
mass. However, the age of an individual system is a notoriously 
difficult parameter to ascertain and so we will return to discuss 



the reliability of this age in Sect. 7.3 



5. Luminosity Determination 



McCaughrean et al. ( 2004 1 used their photometry of s Indi Ba, 
Bb and estimated bolometric corrections from the spectral type- 
Mb i relation of |Golimowski e t al. (2004) to derive luminosities 
of log L/L.0 = -4.71 and -5.35 for s Indi Ba and Bb, respec- 
tively, with estimated uncertainties of ±20%. Here we use our 
flux calibrated spectra to derive luminosities directly following 
a similar process to Golimowski et al. (2004). 

We sum the 0.63-5.1 /an flux with linear interpolation in the 
unobserved regions 2.5-2.86 //m and 4.19-4.53 fim. This is ex- 
tended to the mid-IR using the Spitzer Infrared Spectrograph 
(IRS, |Roellig et al.||1998|) 5-15//m (|Roellig et al.|[2004] i and 



10- 19 /mi spectra (|Mainzer et al. 2007) of the unresolved s Indi 
Ba, Bb system. To approximate the mid-IR spectrum of each 
object and to ensure an accurate representation of the absolute 
fluxes, we also used the resolved mid-IR VLT/VISIR photome- 
try of Sterzik et al. (2005 i. These two determinations of the ab- 



solute flux of the combined system differ significantly. We found 
that the VISIR photometry would have to be scaled upward by 
~50% to be fully consistent with the Spitzer spectrum. We there- 
fore used these two measurements as the bounds on the total flux 
in this region and used the measured flux ratio of ~2.1 at 8.6 pm 
and 10.5 pm from Sterz ik et al.| ( |2005] ) to produce approximate 
spectra for each object for our luminosity determination. 

Beyond this we used a blackbody tail with T e ff given by 
the model fit to the spectra (see Sect.[8]). This however con- 
tributes only 0.5% and 0.9% to the total flux for e Indi Ba and 
Bb, respectively, so the precise temperature used is unimportant. 
The mid-IR spectrum used to extend our spectroscopic measure- 
ments contributes ~8.5% and ~14.5% of the total flux of s Indi 
Ba and Bb, respectively. We did not attempt to extend our ob- 
served spectra blueward for the luminosity derivation as our V- 
band photometry shows the flux to have dropped by 2-3 orders 
of magnitude relative to the /-band flux for both s Indi Ba and 
Bb. 
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The resulting luminosities are log L/L0 = -4.699 + 0.017 
and -5.232 ± 0.020 for s Indi Ba and Bb, respectively, with 
the uncertainties derived from our photometry, 20% uncertainty 
on the mid-IR fluxes, and an assumed 3% uncertainty on the 
absolute flux of the Vega spectrum used for flux calibration 
( |Mountain et al.|1985| |Hayes|1985] ), all of which dominates the 
distance uncertainty. The difference between our determination 
of the luminosity of s Indi Bb and that of McCaughrean et al. 
2004 ) is likely explained by the increased scatter in the spectral 



type-BC/;- relation of Golimowski et al. ( 2004 1 at late-T spectral 
types. 



6. Dynamical Masses 

Ongoing relative astrometric monitoring of the s Indi Ba, 
Bb orbit since May 2004 has allowed a preliminary system 



mass of 121 + lMj U p to be determined (Cardoso et al. 2009a 
McCaughrean et alT||2009| in prep.). As mentioned previously, 



absolute astrometric monitoring is also ongoing (Cardoso et al. 
2009b in prep.) which will determine the mass ratio of the sys- 
tem allowing a model independent determination of the individ- 
ual masses of both s Indi Ba and Bb. However, we can already 
provide some constraints on the individual masses of the brown 
dwarfs on the basis of our photometric and spectroscopic obser- 
vations. Since we know these objects to be physically bound and 
so, in all likelihood, co-eval, e Indi Ba must be more massive 
than its fainter companion. Additionally, we also know from the 
low temperatures required to produce T dwarf spectra that both 
objects must be substellar, therefore the more massive compo- 
nent must have a mass below the hydrogen burning minimum 
mass (HBMM) of -0.070 M© (73 M Jup ) for a cloudy, approxi- 
mately solar metallicity source (Chabrier et al. 2000; Saumon & 
Marley|2008 ' 



Together then, the masses of these two brown dwarfs are con- 
strained to be between 73 + 47 Mj up and 60 + 60 Mj UD , with the 



IJup, 



latter being unlikely due to the different luminosities of the two 
brown dwarfs. Therefore, the mass of s Indi Ba must be in the 
range 60-73 Mj up and s Indi Bb in the range 47-60 Mj up , but 
with a sum of any pairing equal to 121+1 Mj up . 



7. Evolutionary Model Comparisons 

7.1. photometry 

In Table [2] we presented the apparent magnitudes of both com- 
ponents of the s Indi Ba, Bb system in the FORS2 VRIz, and 
MKO JHKL'M' filters. Here we will compare the photometric 
predictions of the solar metallicity COND03 evolutionary mod- 
els ( Baraffe et al.||2003 hereafter B03) with the absolute mag- 
nitudes derived from our observations. However, before doing 
so, it is important to note the limitations of these models even 
if they represent the current state-of-the-art. Although reason- 
able predictions can be made for the near-IR colours of mid- 
late T dwarfs, the colours of early T dwarfs fall somewhere be- 
tween the predictions o f the COND03 and the DUSTY00 models 
( |Chabrier et"ai^|200"0"l >. The COND03 evolutionary models ne- 
glect the effect of dust opacity in the radiative transfer, whereas 
the DUSTY00 models included dust but once formed it remained 
in the atmosphere. The newer BT-Settl atmosphere models ac- 
count for the settling of some species from the atmosphere, so 
once incorporated into the evolutionary models may provide a 
more realistic match to observed colours across the L and T spec- 
tral classes. 
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Fig. 12. M K -(J - K) colour-magnitude diagram with the 5 Gyr 
COND03 isochrone ( jBaraffe et al.|2003[ > with masses (in M ©) 



indicated as filled squares. The DUSTY00 5 Gyr isochrone 
( Chabrier et al.|2000| l is also shown (upper, dashed line). The 



positions of e Indi Ba and Bb are marked with filled circles, the 
filled triangles show the T dwarf observations of |Knapp et al.| 
(2004| and Chiu et al. (2006), and the metal -poor T6 subdwarf 



2MASS 0937+2931 is labelled. All magnitudes are in the MKO 
system. We do not show isochrones for different ages as they 



show little variation in position, except that the same masses cor 


respond to different 


magnitudes. 
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Fig. 13. Mk-(K — L') colour-magnitude diagram with the 5 Gyr 
COND03 isochrone with masses (in Mq) indicated as filled 
squares. The DUSTY00 5 Gyr isochrone is also shown (upper, 
dashed line). Symbols are the same as in Fig. [72] w ith the unla- 
beled data from 



Golimowski et a 



Knap p et al.| ([2004ji, |Chiu et al.| ( [2006| l, and 
(2004|l. The data are in the MKO system, 



whereas the model L' magnitude is in the Johnson-Glass system. 



Figures 12 13 



and 14 show colour-magnitude diagrams 
(CMDs) comparing our s Indi Ba, Bb photometry with the 
COND03 and DUSTY00 5 Gyr isochrones along with other T 
dwarf observations from the literature. The model JHK magni- 
tudes have been transformed from the CIT system to the MKO 
system using the colour relations of Stephens & Leggett (2004), 
while the model L' magnitudes are left in the Johnson-Glass sys- 
tem and the / magnitudes in the Bessell system. 
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Fig. 14. M - J) colour-magnitude diagram with the 5 Gyr 
COND03 isochrone with masses (in Mq) indicated as filled 
squares. The DUSTY00 5 Gyr isochrone is also shown (dashed 
line). Symbols are the same as in Fig. [12] /-magnitudes are in the 
MKO system and /-magnitudes use the FORS2 Bessell / filter. 



In Fig. 12 s Indi Ba is seen to lie ~0. 5-1.0 mag redward 
of the model isochrone, while e Indi Bb is ~0.5 mag blueward. 
The position of the metal-poor (-0.4 < [M/H] < -0.1) T6 subd- 



warf, 2MASS 0937+2931 (Burga sser et al.|2003[) is m arked and 
seen to be bluer than e Indi Bb. |Schilbach et aLf ( |2009[ > find that 
2MASS 0937+2931 falls above an extrapolated [M/H] =-0.5 
IBaraffe et ai] ( fl998] l isochrone in an M K -{J - K) CMD and is 
possibly part of the thick disc or halo population. By contrast, 
the position of s Indi Bb suggests only a slightly sub-solar metal- 
licity. 



Figure 13 shows that the predicted (K - L') colours are a 
better match to our and other early T dwarf observations as the 
effects of clouds are much reduced in this region, but observed 
colours are still consistently bluer than the models implying that 
spectral features in these regions are not well reproduced by 
these models. 

When optical colours are used in a CMD, the mismatch be- 
tween observations and models increases. For late T dwarfs such 



as e Indi Bb, as shown in Fig. 14 our observed (I-J) colours are 
significantly redder than model predictions even after accounting 
for differences due to the use of different filter systems. These is- 
sues will be addressed when we compare our observations with 
the more up-to-date atmospheric models in Sect.JH] 

This discrepancy between the isochrones and the observa- 
tions for the early T dwarfs is an indication of the complexity 
of their atmospheres where the role of clouds is greater than for 
later types. As brown dwarfs transition from dusty L to cloud- 
free mid-late T dwarfs, the observations are found to he between 
the COND03 and DUSTY00 models ( |Baraffe et al.l2003] i. 

7.2. physical properties 

Although the previous section has shown that neither the 
COND03 nor the DUSTY00 atmospheres necessarily reproduce 
the observed colours of T dwarfs faithfully, the comparison of 
properties such as luminosity, mass, and effective temperature 
should be more reliable as they are not so strongly dependent on 
the specifics of the atmospheric model used. That said, Chabrier 



Table 7. Predictions of the parameters of s Indi Ba and Bb from 
the COND03 evolutionary models. These parameters are derived 
for the age range of 3.7-4.3 Gyr given by the measured total 
mass (Ba+Bb) and individual luminosities. 



Source 


Mass 


Radius 


T c ff 


logg 




(M Jup ) 


(R©) 


(K) 


(cms -2 ) 


£ Indi Ba 


67.6-69.1 


0.080-0.081 


1352-1385 


5.43-5.45 


£ Indi Bb 


50.0-54.5 


0.082-0.083 


976-1011 


5.27-5.33 



140 



120 



100 - 




et al. ( 2000 1 showed that the luminosity and effective temper- 



ature evolution at a specific mass are not entirely independent 



Fig. 15. Variation of mass with age for lines of constant lumi- 
nosity interpolated from a fine grid of COND03 evolutionary 
models for the observed luminosities of s Indi Ba and Bb (lower 
two curves). The upper two curves show the predicted combined 
mass at each age accounting for the uncertainties on the individ- 
ual luminosities, while the constraints on the system mass from 
the dynamical monitoring are indicated by the upper two hori- 
zontal dotted lines. The lower horizontal dashed line indicates 
the hydrogen burning minimum mass of 0.070 M0 and the 3.7- 
4.3 Gyr age predicted by the COND03 models is shown by the 
two vertical dotted lines. 



of the treatment of dust in the atmosphere. They find the differ- 
ence between the DUSTY00 and COND03 model predictions, 
for a given age, can be up to 10% in effective temperature and up 
to 25% in luminosity with the COND03 models predicting sys- 
tematically higher effective temperatures than the corresponding 
DUSTY00 models. With the precision of our mass and luminos- 
ity determinations, these uncertainties are significant. However, 
for the age and mass range of the s Indi Ba, Bb system, the 
DUSTY00 evolution models predict only slightly less efficient 
cooling than the COND03 models, resulting in ~3% smaller 
masses, 5% larger radii and less than 3% differences in effec- 
tive temperature than the COND03 models. 

Additionally, Saumon & Marley (2008) have recently com- 
pared their evolutionary model predictions with the Lyon 
DUSTY00 and COND03 models, concluding that the differ- 
ences between the corresponding cloudy and cloud-free models 
at intermediate ages can be explained by the lack of electron 
conduction in the |Saumon & Marley (2008) models. The differ- 
ences seen when compared to the |Burrows et aT ( 1997| ) models is 
found to be mainly due to the use of older non-grey atmospheres 
in the Burrows et al. ( 1997) models. We therefore employ only 
the Lyon models in our comparison. 

Given the observed luminosities of s Indi Ba and Bb, the 
evolutionary models can be used to predict the corresponding 
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Table 8. Physical parameters of s Indi Ba, Bb derived using the observed luminosities and the COND03 models of Baraffe et al. 
p003] l for three ages: 1, 5, and lOGyr. 



Source 


logL/L© 




Mass 
(M Jup ) 






Radius 

(Ro) 






T c ff 
(K) 






(cms 4 ) 






IGyr 


5Gyr 


lOGyr 


IGyr 


5Gyr 


lOGyr 


IGyr 


5Gyr 


lOGyr 


IGyr 


5Gyr 


lOGyr 


e Indi Ba 


-4.699 


42.0 


69.8 


74.1 


0.093 


0.080 


0.079 


1274 


1387 


1393 


5.10 


5.45 


5.49 


e Indi Bb 


-5.232 


28.0 


57.2 


66.2 


0.097 


0.081 


0.077 


929 


1019 


1051 


4.88 


5.36 


5.47 



masses as a function of the system age, as shown in Fig. 15 It 
is immediately apparent that the original age estimate of 0.8- 
2.0 Gyr is inconsistent with the evolutionary models given the 
newly established mass of 121±1 Mj up . From the COND03 mod- 
els, we find that an age range of 3.7-4.3 Gyr is necessary to ac- 
commodate the mass and luminosity constraints. Table [7] shows 
that for this age range, the COND03 models predict effective 
temperatures of 1352-1385 K and 976-1011 K for e Indi Ba 
and Bb, respectively, and a mass ratio in the range 0.73-0.78. 
Similarly, Table[8] shows the predicted physical parameters of e 
Indi Ba, and Bb for the observed luminosities and ages of 1, 5, 
and 10 Gyr. 

7.3. age of the s Indi system 

This age range of 3.7-4.3 Gyr for s Indi Ba, Bb is significantly 



larger than the age of 0.8-2.0 Gyr estimated by |Lachaume et al 
(|1999|l fo r e Indi A based on the rotational period given by |Saar 



|& Osten| ( |1997[ >. In fact, the measured v sin i for s Indi A was too 
small to derive any meaningful rotational period and this period 
was inferred from the Ca II H&K emission observed by |Henry| 
et al. ( 1996), rendering it less reliable. Barnes (2007 ) also calcu- 



lates an age of 1.03 ±0.1 3 Gyr from rotation, although this uses 
the same inferred 22 day period. Lachau me et al.| (fl999) also 
derive an age from the Ca II H&K activity using their R^ K -age 
relation, which suggests an age in the range 1-2.7 Gyr. However, 
in choosing a young age based on this activity indicator, they ig- 
nored the much greater estimate of >7.4 Gyr which they derived 
from the kinematic properties of the system. 

Age constraints are also available if one considers the larger 
moving group of which s Indi A is the eponymous member, and 
for which Cannon ( 1970 1 quoted an age of 5 Gyr derived from 
the observed (B - V) colour of an apparent red giant clump. 
However with only 15 members (Egge np958[ ), it is not clear that 
there would be enough objects at the end of their main-sequence 
lives to fit this reliably. Furthermore, the moving group members 
given in Eggen]( [T958| l were presumably used in Cannon ( 1970} , 
although no details of the members used are given. This mem- 
bership list was later revised by Eggen ( 1971| l which excluded 
five stars from the original list of members. Thus, the age of this 
apparently elderly moving group is not a strong age constraint by 
itself. However A Aurigae, another member of the moving group, 
has an isochronal age of 5.8+0.43 Gyr from Soubira n & Girard] 
(2005) which adds support to an intermediate age for the entire 
moving group, although this technique has been shown to be un- 
reliable when applied to individual objects. This worsens toward 
later spectral types where the isochrones become degenerate, so 
isochrone fitting would be unreliable for e Indi A, but in the case 
of A Aurigae the isochrones are relatively well separated. 

(2002) classifies s Indi A as a 



Finally, Rocha-Pinto et al. 



"chromospherically young, kinematically old" star based on the 
disagreement between the age from activity indicators and the 
observed space velocities. They find a chromospheric age of 
0.39 Gyr, but note that s Indi A has no obvious Li I absorption 



line and so argue that it is an older star than activity suggests. 
The origin of these stars is as yet unknown. 

After further consideration then, an age of ~3.7-4.3 Gyr for s 



Indi Ba, Bb as predicted by the evolutionary models in Sect. 7.2 



seems plausible at least. The discrepancy between the various 
age estimates highlights the problem of deriving ages for indi- 
vidual objects, as any one indicator cannot be entirely trusted. 
To address this issue, we are currently investigating the feasibil- 
ity of obtaining high resolution, time-resolved spectroscopy of s 
Indi A to derive an age via asteroseismology. 



8. Atmospheric Model Comparison 

The comparisons presented in Sect. |7.1| demonstrated that the 
COND03 and DUSTY00 models yield a relatively poor match 
to the observed photometric properties of s Indi Ba and Bb. 



Using the newer BT-Settl atmosphere models (Allard et al. 2003 



Allard 2009 in prep.), which account both for the formation 
and optical effects of dust in the atmosphere and for settling 
of condensates under steady-state conditions, the calculation of 
synthetic spectra and colours can reproduce these observations 
much more accurately. By allowing the description of partially 
settled clouds, these model atmospheres are a more appropriate 
comparison in particular for early T dwarfs, where the atmo- 
spheres transition from being dust-dominated to cloud-free, and 
so here we compare BT-Settl atmosphere models with our de- 
tailed spectral observations of e Indi Ba and Bb. 

These atmosphere models have yet to be incorporated into a 
self-consistent set of interior structure and evolution models, but 
since even the differences between the DUSTY00 and the com- 
pletely cloud-free COND03 models result in only small discrep- 
ancies in cooling rates and radius evolution, our analysis should 
not be much compromised by applying the newer models on top 
of the older evolution models. 

The BT-Settl models are based on version 15 of the PHOENIX 
stellar atmosphere code (Hausch ildt & Baron| [l999) with up- 
dated opacity databases, including among others the recent water 
line list of Barber et al. ( 2006 ) and extended methane line lists 
( |Homeier et al.||2003| l. The chemical equilibrium code used in 
the DUSTY 00 and COND03 models ( |Allard et al.|200Tl|Baraffe| 
et al.||2003 1 has been modified to include grain settling effects 
( |Allard et al.||2003] l. At each layer in the model atmosphere, 
the dust grain number densities are calculated in equilibrium to 
the gas phase. The timescales for condensation and growth pro- 
cesses and gravitational settling are then calculated following 
Rossow ( 1978 ) and compared to the turbulent mixing timescale 
to predict the median size and the fraction of grains which have 
settled. This fraction is then removed from the composition 
and a new equilibrium obtained, iterating the process until the 
grain density no longer changes. For more detailed results and 
a comparison with other cloud models see Helling et al. (2008 ). 
Turbulent mixing in the present models is calculated by interpo- 
lation from 2D and 3D radiation hydrodynamic (RHD) models 
with the C05BOLD code, extending the results of Ludwig et al. 
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Wavelength (|im) 

Fig. 16. Near- to thermal-IR spectrum of e Indi Ba (black line) with sub-solar metallicity ([M/H]=-0.2) BT-Settl spectra with 
T e ff=1300, 1320, 1340K and log g=5.50 (green, red, and blue lines respectively). The T e ff=1240K spectrum (magenta line) is also 
shown to indicate the large difference between the observations and spectral models for this lower effective temperature . All spectra 
have been smoothed to 60 A FWHM and the observed spectrum median filtered to remove the lowest signal-to-noise points. 



(2006J to lower temperatures by including a self-consistent dust 
module (Freytag et al. 2009). As such, the BT-Settl models do 
not include any adjustable parameters, except for some freedom 
in the translation of the hydrodynamic velocity field in the con- 
vective overshoot region into an effective timescale. This conver- 
sion has been chosen such as to optimally reproduce the entire 
photometric sequence from early-L to mid-T dwarfs. 

For each model spectrum of known effective temperature, 
we inferred a radius using our observed luminosities of s Indi 
Ba and Bb, which along with the known distance, allowed us to 
convert the emergent flux density of the atmospheric models to 
the flux that would be observed from Earth. Therefore, for each 
model spectrum we were able to compare predicted absolute flux 
levels against those observed without reliance on evolutionary 
models. In other words, we were not free to normalise the model 
spectra to match our observations. 

The model grid sampled effective temperature in steps of 
20 K, surface gravity in steps of 0.25 dex, for solar and slightly 
sub-solar metallicity ([M/H]= 0.0, -0.2). The metal abundances 
of Grevesse & Noels ( 1993| l were employed instead of the more 
recent Asplund et al. (2005) abundances as the validity of the lat- 
ter determination of oxygen abundance is the subject of ongoing 
debate given its effect on previously well-matched helioseismo- 
logical theory and observations (cf. Ayres|[2008{ |Caffau et al.| 
|20T)8] l. 

8.1. effective temperature effects 



Figure 16 shows the BT-Settl near-IR spectral models with ef- 
fective temperatures in the range 1300-1340K and log g=5.50, 
[M/H]=-0.2 compared to our spectrum of s Indi Ba. Since 



Kasper et al. (2009) find acceptable fits to low resolution near- 
IR spectra of s Indi Ba using models with effective temperatures 
of 1250K and 1300 K, we include our 1240K model to show 
the large mismatch such a low effective temperature would have 
with the BT-Settl models. Each of the models have an absolute 
flux scale set by the effective temperature and the observed lumi- 
nosity of the corresponding object, therefore there is no scaling 
of model spectra to improve the fit to observations. We find ef- 
fective temperatures in the range 1300-1340 K produce the most 
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Fig. 17. 1 .23-1 .27 fim spectrum of s Indi Ba (thicker line, fourth 
from top) and the [M/H]=-0.2, log g=5.50 BT-Settl spectra with 
T e ff=1280-1380K (increasing bottom to top in steps of 20 K). 
All spectra have been normalised and offset for clarity. Model 
spectra have been smoothed to 2.4 A FWHM. 



reasonable fit to the flux level across the spectrum and to indi- 
vidual features. 

For higher effective temperature models, the KI lines at 
1 .25 fim are too strong and the FeH and CrH features around 
1 fim are too deep, while by 1280 K the 1 fim features are no 
longer reproduced and the depth of the KI doublet is signifi- 
cantly reduced (see Fig. 17 i. Additionally, at lower temperatures, 
we find that the flux level of the near-IR peaks become increas- 
ingly difficult to reconcile with the observations. 



Figure 18 shows the BT-Settl near-IR spectral models with 
effective temperatures in the range 880-960 K and log g=5.25, 
[M/H]=-0.2 compared to the data from s Indi Bb. While the op- 
tical spectrum is problematic at any effective temperature (see 
Sect. |8.4| >, there is relatively little variation within this effective 
temperature range in the thermal-IR. Based on the near-IR alone 
then, we find the most reasonable fits to have effective tempera- 
tures in the range 880-940 K. Figure [18] shows that at tempera- 
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Wavelength (^im) 

Fig. 18. Near-IR spectrum of s Indi Bb (black line) with sub- 
solar metallicity ([M/H]=-0.2) BT-Settl spectra with T eff =880, 
900, 920, 940, 960 K and log g=5.25 (cyan, green, red, blue, 
and magenta lines respectively). All spectra have been smoothed 
to 60 A FWHM and the observed spectrum median filtered to 
remove the lowest signal-to-noise points. 
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Fig. 19. 1.23-1. 27 fim spectrum of s Indi Bb (thicker line, sec- 
ond from bottom) and the [M/H]=-0.2, log g=5.25 BT-Settl 
spectra with T e ff =840-960 K (increasing bottom to top in steps 
of 40 K). All spectra have been normalised and offset for clar- 
ity. All spectra have been smoothed to 10 A FWHM for easier 
comparison. 



tures below 900 K, the flux level of the 1.1 fim and the 2.1 fim 
peaks becomes progressively more depressed, and as seen in 
Fig. [19] at higher temperatures, the 1.25 /mi KI lines are con- 
siderably deeper than observed, although the shapes of the tops 
of the 1 . 1 /im and 1 .25 fim peaks are not well-matched by any of 
the models. 



8.2. metallicity and surface gravity effects 

In low-mass stars, the effect of decreasing metallicity is to in- 



crease the effective temperature at constant mass (Baraffe et al 

[1997]! 



with the magnitude of the effect decreasing toward the 
substellar boundary. However, it is not clear what happens in the 
lower-mass regime where degeneracy effects may alter the sit- 
uation. Since no evolutionary models exist for sub-solar metal- 



licity, substellar objects, we cannot provide a quantitative analy- 
sis of the evolutionary effect of slightly sub-solar metallicity. As 
discussed in Sect.[4] studies of the parent main-sequence star, s 
Indi A, derive a metallicity in the range [Fe/H]=-0.23 — 1-0.06, 
with the most recent study supporting [Fe/H]=-0.2. The effects 
of lower metallicity and increasing surface gravity are somewhat 
complementary, so we compare our observations with BT-Settl 
models of surface gravity from log g=5.00 to 5.50 and with solar 
([M/H]=0.0) and slightly sub-solar metallicity ([M/H]=-0.2). 

Figure 



20 



shows the near- to thermal-IR spectrum of s 
Indi Ba compared to models with an effective temperature of 
1320 K, surface gravities log g=5.25 and 5.50, and metallici- 
ties of [M/H]=0.0 and -0.2. The effect of higher metallicity 
and lower surface gravity is shown to be most prominent in the 
J- and //-bands. The flux in the 1.1, 1.25, and 1.6 /an peaks 
is suppressed, while the 2.2/on peak flux is over-estimated. 
Additionally, the shape of the K- and L-bands is inconsistent 
with the observed spectrum. We find the higher surface gravity 
(log g=5.50) and sub-solar metallicity ([M/H]=-0.2) to be the 
better fit to the spectrum of s Indi Ba. 



Figure 21 shows the near-IR spectrum of s Indi Bb compared 
to models with an effective temperature of 920 K, surface gravi- 
ties log g=5.00 and 5.25 cm s~ 2 , and metallicities of [M/H]=0.0 
and -0.2. We neglect the longer wavelength data here as the 
models show little variation. As with the comparison of models 
with different effective temperatures, we find the //-band is rela- 
tively invariant to the different surface gravities and metallicity. 
The effects of surface gravity and metallicity are most apparent 
in the J- and Zf -bands. The general trend for dust-free T dwarfs 
such as s Indi Bb is well established (cf. Legg ett et al.| |2009): 
both the lower metallicity and higher gravity result in higher at- 
mospheric pressures, favouring the formation of methane and 
increasing collision induced absorption (CIA) opacity around 
2 fim. 

The solar metallicity model with surface gravity log g=5.00 
is the most inconsistent with our observations as the 1.1 
and 1 .25 /im peaks are significantly under-estimated while the 
2.1 /an peak is significantly over-estimated. As explained ear- 
lier (see Sect. [8), we have chosen to use the Grevesse & Noels| 
(1993) abundances which is likely the cause of the differences 
seen between our solar metallicity model spectra and those used 
in Kasper et al. (2009]). The other models make very simi- 
lar predictions to one another, except in the ZT-band. Although 
from Fig. 21 the log g=5.00, [M/H]=-0.2 model is marginally 



favoured, the comparison to all models within the effective tem- 
perature range of 880-940 K favours log g=5.25, [M/H]=-0.2. 
Additionally, none of the models can reproduce the lower flux 
seen in the 2.2-2.5 /on region. 

In summary, we find BT-Settl atmosphere models with ef- 
fective temperatures in the range 1300-1340 K, surface gravity 
of log g=5.50, and slightly sub-solar metallicity of [M/H]=-0.2 
are the most reasonable fits to the observed optical to thermal-IR 



spectrum of s Indi Ba. In Fig. 22 we show the 1300 K spectral 



model compared to the data. It is clear that this model does not 
provide a perfect match to the near-IR peaks, however general 
features are reproduced. The inset plot shows the mismatch in 
the shape of the K I profile at 0.77 /im which may suppress the 
continuum to beyond 1 fim (Burr ows et al.|[2000) ) and so affect 
the relative strengths of the 1 . 1 /an and 1 .25 fim peaks. Although 
the flux in the 2.8-3.2 fim region is underestimated, the 3.3- 
4.2 fim region is well-reproduced. However, in the M-band the 
model predicts higher fluxes than observed, suggesting insuffi- 
cient understanding of the sources of opacity. 
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Wavelength (pin) 

Fig. 22. Optical to thermal-IR spectrum of s Indi Ba (black line) and the best-fit BT-Settl spectrum (red line, T e ff=1300K, log 
g=5.50, [M/H]=-0.2). Here the optical and near-IR data have been smoothed to 17 A FWHM and median filtered to remove the 
lowest signal-to-noise points. The model has been smoothed to 60 A FWHM beyond 2.5 fim. Inset: the optical portion of the spectra 
smoothed to 6.5 A FWHM shown with the flux on a logarithmic scale. Also shown is a model fit to the optical spectrum with depleted 
alkali abundances (blue line). For this early T-dwarf, the standard alkali abundances provide the better match to the observed flux 
levels and spectral morphology. Additionally, the lithium abundance has been reduced by a factor of 1000 to allow a comparison 
with the observations which show no detectable lithium at 6707 A. 
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Fig. 23. Optical to thermal-IR spectrum of s Indi Bb (black line) and the best-fit BT-Settl spectrum (red line, T e ff =920 K, log g=5.25, 
[M/H]=-0.2). The spectra have been smoothed as in Fig. 22 For this late T-dwarf, the depleted alkali model provides the better 
match to the observed flux levels and spectral morphology in the optical which may be an indication of the formation of feldspars 
in late-T dwarfs which is not evident in earlier types. No lithium depletion has been implemented in these models. 



For s Indi Bb, we find the BT-Settl models with effective 
temperatures of 880-940 K, surface gravity of log g=5.25, and 
metallicity of [M/H]=-0.2 provide the best match to the obser- 
vations. Figure [23] shows the 920 K model does not provide a 
perfect match to the near-IR peaks, although this is less pro- 
nounced than for Ba. The inset plot shows the large mismatch in 
the shape and flux level of the K I profile at 0.77 fim and the ef- 
fect of employing a model with depleted alkali absorption. This 
effect is also evident in the optical colour-magnitude diagram of 
Fig. [14] where s Indi Bb falls -1.5 mag redward of the COND03 



evolutionary models. We return to discuss the cause of this large 
effect in Sect. 8.4 The mismatch around 1.62-1. 74 fim is mainly 



due to incomplete knowledge of CH4 absorption in this range. 
In the 2.8-3.5 fim region, the shape of the spectrum due to CH4 
absorption is not well fit. In particular, the drop in flux predicted 



at ~4. 1 fim is not observed, and the M-band flux is again over- 
estimated. 

The KI-H2 quasi-molecular satellite feature at ~0.69 fim dis- 
cussed in Alla rd etliL] ( |2"00 7) is seen here in both s Indi Ba and 
Bb (inset plots of Figs. [22] and [23 1. The shape of the feature is 
not perfectly reproduced, mainly because the effects of KI-He 
absorption have yet to be included. To a lesser extent, the ill- 
fitting wings of the very wide, pressure-broadened absorption 
lines of K I at 0.77 fim and Na I at 0.59 fim may also play a role. 



8.3. unidentified feature at 1.35-1. 40 fim 

We note a spectral feature at ~ 1.35-1 .40 pm in the spectra of 
both e Indi Ba and Bb which is poorly reproduced by the BT- 
Settl models (see Figs. [22] and [23]) and has not been identified 
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Fig. 20. Near- to thermal-IR spectra of s Indi Ba (black line), the 
BT-Settl spectrum (T eff =1320K, log g=5.50, [M/H]=-0.2, red 
line), the same with log g=5.25 (blue line), and with log g=5.50 
but with solar metallicity ([M/H]=0.0, green line). All spectra 
have been smoothed to 60 A FWHM and the observed spectrum 
median filtered to remove the lowest signal-to-noise points. 
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Fig. 21. Near-IR spectra of s Indi Bb (black line), the BT-Settl 
spectrum (T eff =920K, log g=5.25, [M/H]=-0.2, red line), the 
same with log g=5.00 (blue line), with log g=5.00, [M/H]=0.0 
(magenta line), and log g=5.25, [M/H]=0.0 and (green line). All 
spectra have been smoothed to 60 A FWHM and the observed 
spectrum median filtered to remove the lowest signal-to-noise 
points. The thermal-IR spectra have not been shown as there is 
little to distinguish between models. 



in all spectra of the T dwarf standard stars. Many observers 
remove this region from published spectra due to the high tel- 
luric absorption and so there are few available comparisons. 
Nevertheless, the region contains valuable information on the 
continuum flux level in these deep absorption bands present in 
T dwarfs. A similar feature is however seen in the spectrum 
of the Tl spectral standard SDSS0151 + 1244 ( |Burgasser eTaL] 
2006 b]), the T8.5 and T9 dwarfs ULAS1238 and ULAS1335 
(Burningham et al. 2008 1, and some L dwarfs (e.g. 2MASS 
J1507-1627(L5), |Burgasser|2007| l. 

While this feature could be caused be problems with the tel- 
luric correction in this region, we believe that it may be intrinsic 
to these objects. We are satisfied that this feature is not an arte- 



fact of our spectral extraction. The rise in flux occurs part way 
through the last of the seven spectral orders which were com- 
bined to produce our /-band spectrum (see Appendix |C|), so the 
shape is not due to a scaling mismatch between adjoining re- 
gions. This region corresponds to the gaps in our standard star 
spectrum where high telluric absorption required replacement by 
the solar model of Kurucz (see Sect. |2.6[ ), but this is not the cause 
of the feature, as the solar model has only weak spectral features 
here. 

The BT-Settl model for s Indi Ba shows some structure in 
the deep water bands between the J and H peaks that is qual- 
itatively similar to the observed feature, though the flux is un- 
derestimated by a factor of ~2. In the s Indi Bb model the dis- 
agreement is much worse, with about an order of magnitude mis- 
match between the modelled feature (which is nearly invisible 
at the scale of Fig. 23 1 and the observation. Still, this suggests 
that the feature is due to the structure of the strongest part of 
the water absorption bands, which in this part of the spectrum 
form fairly high up in the atmosphere, corresponding to opti- 
cal depths of Tp> osse i an d=0. 1-1 .0 and temperatures of 700-1000 K. 
Incidentally, two other features which also have systematically 
underestimated flux in the present models form at a similar level, 
namely the collisionally-induced absorption (CIA) capping the 
flux peak in the K band, and the wings of the potassium dou- 



blet centred on 0.77 /mi (see Sect. 8.4 1. A possible explanation 



for the mismatch in all these cases might be an underestimated 
local temperature, and thus source function, at this atmospheric 
level. A toy model indicates that an increase in temperature of 
200^400 K could reconcile the modelled and observed flux lev- 
els, but at this point we have no reasonable idea for the cause 
of such heating (e.g. for an additional opacity source causing a 
corresponding back-warming). A full explanation of this feature 
must also account for the apparent difference in strength between 
objects of the same spectral type. 

8.4. alkali depletion 

The optical spectra of T dwarfs are dominated by alkali reso- 
nance lines which are the result of a balance between the in- 
creased transparency of the atmospheres, due to the progressive 
sedimentation of condensates, and the depletion of alkali metals 
from the gas phase, also due to condensation processes. It has 
been noted ( |Lodders & Fegley|2006) |Burrows|20"09"l > that alkali 
metal depletion cannot occur above T~ 1400 K, below the con- 
densation temperatures of most refractory species and further- 
more, the first condensates are the feldspars ([Na,K]AlSi30g) 
which require aluminium and silicon to form. However, in a 
stratified atmosphere, the latter elements can be expected to have 
already been depleted by higher-temperature condensates before 
feldspars would have a chance to form. In that case, the alkali 
elements could only condense into sulfides and halides (mostly 
Na 2 S, KC1) at temperatures around 1000 K. 

In any case, for mid- to late-type T dwarf atmospheres, the 
depletion of the sodium and potassium is not adequately ac- 
counted for in current models. At the higher temperatures of 
early T dwarfs such as e Indi Ba, the standard BT-Settl model 
reproduces the shape and flux level of the K I doublet at 0.77 /mi 
(and the red edge of the NaD line at 0.589 /mi) reasonably well 
as shown in Fig. [22] (inset), but over-estimates the absorption in 
the later type s Indi Bb by an order of magnitude (Fig. 23 inset, 
red line). 

The line wings of the K I resonance doublet at 0.77 /mi ex- 
tend several 1000 A and thus potentially suppress the flux out 
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Fig. 24. The 6650-6750 A spectra of e Indi Ba and Bb (upper 
and lower thick black lines, respectively) and the BT-Settl mod- 
els with T eff =1300K, log g=5.50, and [M/H]=-0.2 (upper dot- 
ted line) and T eff =920K, log g=5.25, and [M/H]=-0.2 (lower 
dotted line). Both models have been smoothed to 6.5 A FWHM 
to match the observations. The lithium abundance has been re- 
duced by a factor 1000 in the T e ff=1300K model in order to 
provide a meaningful upper limit to the observed levels of Li I 
absorption in the spectrum of s Indi Ba. 



to the 1.1 and even 1.25jum peaks. As pointed out already by 
Burrows et al. ( 2000 1, these line profiles show strong deviations 
from classical Lorentzian wings and thus require a more sophis- 
ticated broadening theory. Detailed spectral models of the alkali 
lines have been successfully used to model the optical spectrum 
of s Indi Ba by |Allard et al.| ( |2007| l, based on improved inter- 
action potentials with H2 and He for the line profiles, and an 
earlier version of the settling framework to calculate the depth- 
dependent abundances of neutral alkali atoms. They were also 
able to identify the quasi-molecular satellite of KI seen in our 
resolved spectra (see Figs. 22 and 23 1. While the overall agree- 



ment of the modelled optical spectrum for s Indi Ba is fair, the 
Bb model, though showing the same general features, underesti- 
mates the observed flux by up to an order of magnitude. Since the 
lines in both brown dwarfs should form at quite similar temper- 
ature levels, it is extremely unlikely that the line profiles are off 
by such a large amount for the cooler component only. One pos- 
sible explanation for the mismatch is that the K I in the gas phase 
is much less abundant than expected from the settling model. 

If we therefore consider additional condensation of the al- 
kali metals into feldspars, which would result in alkali deple- 
tion from the gas phase at somewhat higher temperatures than in 
current models, we can bring the s Indi Bb spectrum into rea- 
sonable agreement with observations as shown in Fig. [23] (inset, 
blue line). However in this case the fit to the s Indi Ba spectrum 
becomes worse. The formation of feldspars implies that Al and 
Si would still have to be present in sufficient quantities at the 
feldspar condensation level and, indeed, recent calculations us- 
ing results from updated RHD simulations (Freytag et al. 12009 



|Homeier et al.|20G9] in prep.) imply that the upmixing of these 
species might be more efficient than previously assumed. Thus 
we suggest that feldspar formation can efficiently deplete these 
species at the temperatures of mid-late T dwarfs. 



8.5. lithium 

The presence and state of lithium in the atmospheres of brown 
dwarfs is governed by mass, present effective temperature, and 
age. The m odels of |D'Antona & Mazzitelli| | |1994) , |Chabrier 



|et al.| ( |1996] l, and |Burke et al. (2004) predict that a brown dwarf 
must have a mass greater than 0.060-0.065 Mp, to be capable of 
reaching the core temperature of ~3xl0 6 K (Bildsten et al. 1997 1 
required for lithium depletion ( 7 Li + p — > 4 He + 4 He). Therefore, 
since objects less massive than ~0.4Mq are fully convective, all 
lithium should have been processed for any brown dwarf above 
0.065 M©. |Chabrier et al. ([1996} predict that an object at this 
mass boundary will deplete its primordial lithium by a factor of 
100 in ~1 Gyr, with faster depletion in higher-mass objects. 

Furthermore, the substellar chemistry models of Lodd ers &| 
Fegley] ([2006} show that Li I is the dominant form of lithium 
down to ~1520K at 1 bar pressure (with the temperature limit 
increasing with increasing pressure), beyond which lithium is 
bound into molecules (LiCl, LiOH, etc.). Kirkpatrick et al 
(2000) and Kirkpatrick et al. (2008) presented spectra of a num 



ber of L dwarfs with and without Li I absorption. They showed 
that the strength of the Li I 6708 A line peaks at L6 and is ob- 
served in some L7 and L8 dwarfs with the strength of the line 
and number of detections decreasing toward later types. 

Here we have higher resolution and higher signal-to-noise 
data and, from the observed spectral type-equivalent width re- 
lation of Kirk patrick et al.| (|2000), we would have expected to 
detect Li I at 6708 A in s Indi Ba if it were present, but not nec- 
essarily in e Indi Bb. Indeed, Burrows et al. (2002) state that 
there is no apparent reason why Li I would not be seen in objects 
as late as T6 with high enough quality data. 

On the other hand, the presence of Li I absorption in our cur- 
rent atmospheric models may be due to unrealistic modelling of 
the removal of solid species from the atmosphere after forma- 
tion (Marley, pers. comm.). As can be seen in the Fig. 24 we 



find no evidence for absorption by monatomic lithium at 6708 A 
in the spectrum of s Indi Ba, while the spectrum of s Indi Bb 
may have inadequate signal-to-noise at this wavelength to judge 
its presence or absence at the low abundance levels expected for 
T6 dwarfs. To approximately quantify the level of lithium deple- 
tion, our models for s Indi Ba must have the proto-solar lithium 
abundance reduced by a factor of at least 1000 (see Fig. 24 1 in 
order to reproduce the data. 

The level of lithium depletion in s Indi Ba is compatible 
with a mass in excess of 0.065 Mq (68 Mj up ). The less lumi- 
nous s Indi Bb on the other hand, must be less massive than this 
limit, and given its lower temperature, has most probably lost its 
atomic lithium to molecules. 



8.6. chemical equilibrium departures 

At the effective temperatures of s Indi Ba and Bb, carbon 
would be expected to be predominantly locked in carbon monox- 
ide (CO) in the hotter, lower layers of the atmosphere, and in 
methane (CH4) in the cooler, upper regions under chemical equi- 
librium (CE) conditions. However, since the detection of the CO 



fundamental band in Gl 229B by|Noll et al. ( 1997), observational 



evidence has accumulated that CO persists in the upper atmo- 
spheres of T dwarfs in excess of its CE abundance. 

|Griffith & YeIIe1 ( [T999| l and |Saumon eTaT[ ( |2003l have shown 
that this excess can be explained by the upmixing of CO from 
the warm deeper layers, since if one assumes sufficiently effi- 
cient turbulent mixing in the upper atmosphere, the kinetic rates 
for the conversion reactions from CO to CH4 are too slow to 
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Fig. 25. The 2.27-2.35 /-im spectrum of s Indi Ba (bottom line) 
showing the CO 2-0 overtone along with models employing dif- 
ferent chemical equilibrium (CE) timescales. All models have 
T eff =1320K, log g=5.50, and [M/H]= -0.2. A chemical equilib- 
rium model is shown at the top, the Prinn & Barshay ( 1977|l 
model second from bottom, and a model following the Yung 



|et al.| ( |1988| l chemistry second from top. All models have been 
smoothed to 4.9 A FWHM to match the observations. The s Indi 
Ba spectrum shows the absolute flux scale, while the Prinn & 
|Barsnay| ( |1977[ ), |Yung et al.] ( |1988[ ) and CE models have been 
offset for clarity by 0.1, 0.8, and 1.0, respectively. 



adjust the mixing ratios to CE abundances. The BT-Settl mod- 
els employ a similar calculation of these CE departure effects, 
but using diffusion coefficients derived from the C05BOLD ra- 
diation hydrodynamics simulations (Freytag et al. 2009 1 rather 



than adjusting them as a free parameter ( |Homeier et al. 2009 
in prep.). The diffusion coefficients thus calculated are there- 
fore height-dependent and not directly comparable to a single 
choice of an eddy diffusion coefficient, with typical values in 
our models ranging between 10 5 and 10 9 cm 2 /s. Since the CE 
departure is most sensitive to the mixing at the temperature level 
relevant for the transition from CO to CH4, i.e. between 1000 
and 1500 K, the diffusion coefficient in this part of the overshoot 
region would best characterise our model, corresponding to 10 7 - 
10 8 cm 2 /s in the s Indi Ba model. In addition, we consider sev- 



eral reaction pathways and timescales besides the one from Prinn 
& Barshay]( |1977[ ), on which the Saumon et al. (2003 ) models are 



based, notably the revised time scale for the scheme of |Prinn & 
Barshay I ( |1977[ l suggested by 



reaction scheme of |Yung et al. 



( 1988) (see also Griffith & Yelle 



Griffith & Yelle (199"9l, and the 



1999)). Among these, the |Yung et al. ( 1988| l model generally 
predicts the fastest conversion rates from CO to CH4, and the 



Prinn & Barshay ( 1977 ) model with the modified rates of Griffith 
& Yelle| ( |1999[ ) the slowest rates, implying the strongest CE de- 



parture effects. 

Since the velocity field derived from the C05BOLD simula- 
tions is used for the description of both the cloud dynamics and 
the CE departures, these models employ a high degree of self- 
consistency, with any change in the mixing properties immedi- 
ately affecting both the dust content and the gas-phase chemistry 
of the atmosphere. On the other hand, other potential uncertain- 
ties in our cloud model would also feed back into the thermal 
profile of the atmosphere and might thus affect the domains of 
different carbon chemistry discussed above. However, any major 
changes in cloud opacity would also inevitably change the spec- 



tral energy distribution and thus produce inconsistencies with the 
observed IR photometry. As the overall fit to the spectral energy 
distribution is good, we thus feel confident that our chemistry 
model is not affected by major uncertainties in the thermal struc- 
ture due to backwarming from the cloud deck. 

As an L-T transition object, marking also the transition from 
CO-dominated to CH4 -dominated chemistry, s Indi Ba is partic- 
ularly sensitive to the reaction details discussed above. The high 
quality of the present K- and L-band spectra thus allow quan- 
titative estimates of the CO and CH4 mixing ratios, enabling 
us to directly test these models of the non-equilibrium chem- 



istry and constrain the relevant timescales. Fig. 25 shows the ob- 
served 2.27-2.35 jum spectrum, where the CO 2-0 overtone band 
(starting at 2.2935 fim) can be identified even on top of the octad 
band of CH4. For comparison, a chemical equilibrium model and 
non-CE model spectra, based on the reaction models of P rinn &| 
Bars hay] (|1977[) and|Yung et al.| (p~988) respectively, are shown. 
The Prinn & Barshay| ( |1977| l model better reproduces the rota- 
tional series of the CO 2-0 R branch extending redwards from 

1988| l models match the 



2.2935 //m, but the CE and Yung et al. 



overall morphology better, including the 2.3 15 jim CH4 band- 
head. However, as seen in Fig. 26 the Prinn & Barshay ( |1977| l 



model does not match the shape of the L-band spectrum, whereas 
the chemical equilibrium model and the Yung et al. ([1988 ) model 
provide reasonable matches to the shape and absolute flux lev- 
els - the CH4 absorption predicted by both models is identical 
within the internal uncertainties of the model atmospheres. Also 



shown is a model employing the |Prinn & Barshay 
with the revised reaction rates of Gri ffith & Yelle 



( 1977 1 model 



1999i which 



shows an even poorer match to the spectral morphology, since in 
this model even less CH4 has formed. However, the depth of the 
CH4 absorption at ~3.3 fim falls between that seen in the spectra 
predicted by the |Yung et"al~| (1988 ) and |Prinn & Barshay | ( fl 977 1 
models, just as the 2.3 fim spectrum indicates a CO abundance 
intermediate between those models. This would suggest that the 
correct reaction timescale might have a value slightly below that 
of Yung et al. ( 1988), or alternatively, the vertical mixing could 



be more efficient than assumed in the BT-Settl models. [Freytag 
|et al. ( 2009) 1 suggest such additional mixing can be produced by 



convectively driven gravity waves, though the mixing efficiency 
of such waves in terms of an equivalent diffusion coefficient is 
still under investigation. 

Our observations provide only limited coverage and resolu- 
tion of the strongest CO absorption feature, the 1-0 fundamen- 
tal band at 4.55 jt/m (see Fig. [7]). However, the resolution of our 
M-band spectra does not allow us to favour any model over the 
other. Although the models differ in the predicted absolute flux 
levels, none of these match the observed flux levels. 

We have shown that the observed Zf-band spectrum of e Indi 
Ba shows evidence of non-CE processes and the L-band spec- 
trum supports a CO-CH4 reaction rate intermediate between that 
predicted by the |Yung et"aT[ (l988) and |Prinn & Barshay|fl977| l 



models. However, we postpone a more detailed analysis of the 
non-chemical equilibrium signatures to a future paper. 

9. Mass limits on further system members 

Using the previously derived flux limits from the deep imaging 
of the system (see Sect. 2.4 1 and an approximate system age of 



5 Gyr, we can derive mass limits of further members using mod- 
els. The COND03 models rule out any system members of mass 
greater than 15M Jup (H-I9.l m ) in the field around e Indi Ba, 
Bb (7-39.7 AU), and also rule out objects more massive than 
34M Jup (ff=15.8 m ) from 7 AU down to 0.8 AU from either Ba 
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Fig. 26. The L-band spectrum of s Indi Ba (black line) and 
different chemistry models with T e ff=1320K, log g=5.50, 
[M/H]=-0.2. A chemical equilibrium model is shown in green, 
the Prinn & Barshay ( 1977 1 model in blue, the Prinn & Barshay 
(|1977|l model with the revised reaction rates of Griffith & Yelle 



( 1999 1 in magenta, and a model following the |Yung et al. 
(1988) chemistry in red. All models have been smoothed to 



60 A FWHM to match observations. 

or Bb. If the system were as young as 1 Gyr, the corresponding 
mass limits would be 6.0 and 13 Mj up , respectively. 

At separations less than 1-2 FWHM (-0.4-0.8 AU), our 
ability to distinguish companions of Ba or Bb is limited by the 
high object flux and pixelation of the PSF. However, if either s 
Indi Ba or Bb had substantial unresolved companions, then their 
combined observed spectra would be notably different compared 
to the spectral standards. The combined spectrum of s Indi Ba 
and Bb is drastically different to that of Ba alone, so any addi- 
tional companion must be significantly fainter, and hence less 
massive, than Bb. 

10. Comparison of Model Predictions 

10.1. previous determinations 

Observations of s Indi Ba, Bb have previously been compared 
to different atmospheric and evolutionary models. |Smith et al.| 
( |2003| ) used R~50 000 spectra in the ranges 1 .553-1 .559 fim and 
2.308-2.3 17 fim to derive an effective temperature for s Indi Ba 
from comparisons of the observed spectra with spectral models 
to fit bands of CO and H2O. They derived effective tempera- 
tures of 1400 K and 1600 K from the two regions and so adopted 
1500 K as the effective temperature of s Indi Ba. As they note, 
this is significantly in excess of that derived in other analyses, 
and the temperatures derived from the two regions are incon- 
sistent. Moreover, the effective temperature of 1600 K derived 
for the 2.308-2.3 17/im region may be affected by the CO/CH4 
chemistry. 



trum which agreed well with the observed combined spectrum. 
However, Sterzik et al. (2005 ) presented mid-IR photometry of 
the individual sources from VLT/VISIR which they suggest is 
not compatible with the absolute fluxes of the individual models 
from Roellig et al. ( |2004[ ). Indeed, as pointed out in Sect.[5] the 
Spitzer and VISIR fluxes differ significantly. Using an assumed 
age of 1 Gyr, |Sterzik et al.| ([2005 ) found an effective tempera- 
ture of HOOK for e Indi Ba. However, our comparison of the 
optical to thermal-IR spectrum excludes such low effective tem- 
peratures. 

Finally, |Kasper et al.| ( |20091) used low res olution near-IR 
spectroscopy and the models of |Burrows et al.| ( |2006[ ) to yield 
effective temperatures of 1250-1 300 K and 875-925 K and sur- 
face gravities of 5.2-5.3 and 4.9-5.1, for e Indi Ba and Bb, re- 
spectively. These temperatures are broadly comparable to ours, 
although in more detail we find the BT-Settl spectral models can- 
not match the observed e Indi Ba spectra at effective tempera- 
tures as low as 1250 K. Their grid of surface gravity was finer 
than used here and they derive values significantly lower than 
ours. For e Indi Ba, we have compared our observations to mod- 
els with log g=5.25 and 5.50, and specifically prefer the higher 
value. Similarly, for s Indi Bb, we have tested models with log 
g=5.00 and 5.25, and again prefer the higher value. The various 
determinations of the effective temperatures of e Indi Ba and Bb 
are summarised in Table|9] 

The situation is complicated by complementary effects on 
the spectral morphology due to surface gravity, metallicity, and 
elemental abundances. We do not believe that the surface gravity 
can be constrained to better than ~0.25 dex. More precise deter- 
minations of the surface gravity can lead to substantially inac- 
curate predictions of the mass as in|Kasper et al.|([2"009[), which 



Burrows 



may do the evolutionary models an injustice. Indeed 
et al. (2006) suggest that the lack of a detailed understanding 
of brown dwarf meteorology may lead to ambiguity in derived 
effective temperatures of ~50-100K and surface gravities and 

-0.3. 

Kasp er et al.| ([2009 ) used their fitted effective temperatures 
and surface gravities to derive ages and masses from the evo- 
lutionary models of |Burrows et al. ( 1997 1. For e Indi Ba, they 
found an age of 1.0-2.3 Gyr and mass of 46-62 Mj llp , and for 
e Indi Bb an age of 1.0-2.0 Gyr and a mass of 29-39 Mj up . 
Although the objects appear to be co-eval and consistent with 
the age estimate of Lachaume et al. (1999), this age range is 
clearly lower than we have derived in the present paper us- 
ing the observed luminosities and dynamically derived system 
mass. Additionally, their predicted total maximum system mass 
(101 Mj up ) is considerably lower than the measured system mass 
(121 Mj up ). This suggests that either there are large systematic 
errors in the evolutionary model predictions, or the effective tem- 
peratures and surface gravities derived from observed spectra are 
inaccurate. This comparison highlights the problems associated 
with deriving the physical properties of other field objects for 
which we do not have a comprehensive set of observations or 
external constraints from a companion star. 



Roellig et al. | <|2004|) and [ Mainzer et al. | (|2007f observed the 1 Q 2 Qur modg/ predMons 



combined spectrum of the e Indi Ba, Bb system using IRS on 
Spitzer. Using the luminosities of |Mc Caughrea n et IT] ( |2004| ) 
and an age of ~ 1 Gyr, they derived evolutionary model param- 
eters of T e ff=1210K, log g=5.10 (for a cloudy evolutionary 
model) and T eff =840K, log g=4.89 for e Indi Ba and Bb, re- 
spectively. Using these parameters along with radii of 0.094 and 
0. 10 R for e Indi Ba and Bb, respectively, |Roellig et aT^pXH] ) 
and |Mainzer et al.| ( |2007| ) generated a composite model spec- 



Our comparison of the observed luminosities and measured total 
mass of s Indi Ba, Bb has allowed us to derive a predicted age 
range of 3.7^1.3 Gyr from the COND03 evolutionary models. 
While this age range is higher than the previously used age esti- 
mate of 0.8-2.0 Gyr from Lachaume et al. (|1999[), it is younger 



than other age indicators for e Indi A in the literature discussed 
in Sect.|7.3| From this age range and the individual luminosities 
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Table 9. Comparison of our model predictions with previous 
studies. The last two entries show the determinations of the ef- 
fective temperatures of the two T dwarfs from this work, both 
from fitting atmospheric models to the observed spectra, and by 
comparing the derived luminosities and dynamical system mass 
with evolutionary models. 



Study 



e Indi Ba T a 



£ Indi Bb T cl 



Smith etal. (2003) 
Roellig et at[|2U04l 
Sterzik et al.lo 
Mainzer et al. 



Kasper et at. (2009) 



This work (atm. models) 
This work (evo. models) 



1400,1600 

1250 

1100 

1250 
1250-1300 
1300-1340 
1352-1385 



800 

800 
875-925 
880-940 
976-1011 



we have extracted the predicted individual masses, effective tem- 
peratures, radii, and surface gravities (see Table|7]). The effective 
temperatures predicted by these models are 1352-1385 K for s 
Indi Ba, and 976-101 1 K for Bb with surface gravities of 5.43- 
5.45 and 5.27-5.33. 

We have also independently derived the effective tempera- 
ture and surface gravity of both objects through a direct com- 
parison of the observed optical to thermal-IR spectra with the 
BT-Settl atmospheric models. For s Indi Ba, we found models 
with effective temperature in the range 1300-1 340 K and sur- 
face gravity of log g=5.50 are the most appropriate, while for s 
Indi Bb we find effective temperatures in the range 880-940 K 
with surface gravity of log g=5.25, both with slightly sub-solar 
metallicity of [M/H]=-0.2. 

We cannot presently reconcile the effective temperatures de- 
rived from atmospheric modelling with those derived from evo- 
lutionary models using the measured system mass at the ob- 
served luminosities. However, it must be noted that the evolu- 
tionary models do not presently incorporate the new BT-Settl at- 
mospheres. The difference between the predicted effective tem- 
peratures for s Indi Bb show the atmospheric models to be 
inconsistent with the evolutionary models. The upper limit on 
the effective temperature from the comparison with atmospheric 
models is 940 K, while the evolutionary models predict effective 
temperatures of at least 975 K. For s Indi Ba, the difference is 
less severe. There is only a 10 K difference between the limits 
on the effective temperature from the evolutionary and atmo- 
spheric models. This is midway between the effective temper- 
ature steps in our grid, although it is clear that by 1360K the 
atmospheric models are no longer consistent with the spectro- 
scopic observations. These differences may be resolved by the 
inclusion of the BT-Settl atmospheric models in the next gener- 
ation of evolutionary models, as the effect on the effective tem- 
perature would presumably be intermediate between that of the 
COND03 an d DUSTY00 models. However, we note that |Dupuy 
|et al.| ( [2009| ) also find that the effective temperatures predicted 
using atmospheric and evolutionary models are in disagreement 
for the components of the L4+L4 binary HD 130948BC. 

The derived surface gravities are consistent between the evo- 
lutionary and atmospheric models. However, due to the grid 
step chosen (0.25 dex) and the complementarity with metallic- 
ity, the surface gravities derived from the comparison of our ob- 
served spectra with the atmospheric models do not provide a 
strong test of the atmospheric and evolutionary model predic- 
tions. Additionally, any mass estimated from the radius deter- 
mined for each spectral model, the fitted surface gravity and our 
derived luminosities, will give a large range of possible masses. 



As described earlier, once the individual masses have been 
determined from the ongoing absolute astrometric monitoring 
|Cardoso et al. ( 2009b[ ), we will be able to test the evolutionary 
models using the individual masses, luminosities, and the same 
age for both T dwarfs. Additionally, once more reliable age de- 
terminations become available we will be able to directly test 
the evolutionary models and determine if the luminosities are 
overestimated for intermediate age brown dwarfs as suggested 
by |Dupuy et"al^ ( |2"00"9"l > for the young system HD 130948BC. 

For the observed luminosity of s Indi Ba and effective tem- 
perature range of 1300-1340 K derived from the comparison to 
atmospheric models, the COND03 evolutionary model predicts 
a mass for s Indi Ba in the range 46-64 Mj llp , and for s Indi Bb 
in the range 16-37 Mj up from its observed luminosity and the ef- 
fective temperature range of 880-940 K. Given the preliminary 
dynamical system mass of 121 ±1 Mj up (Cardoso et al.|2009a|), it 
therefore appears that with current theoretical models and spec- 
troscopically derived effective temperatures, one cannot obtain 
reliable mass predictions for T dwarfs such as these even when 
precise luminosity constraints are available. 

11. Conclusions 

We have presented the results of a comprehensive photomet- 
ric and spectroscopic study of the individual components of the 
nearest known binary brown dwarf system, s Indi Ba, Bb. The 
relative proximity of these Tl and T6 dwarfs to the Earth re- 
sulted in very high quality data, while archival results for the 
well-studied parent star, s Indi A, provide invaluable additional 
information. We find the spectra of these brown dwarfs are best 
matched by the BT-Settl spectral models with T eff = 1300-1 340 K 
and log g=5.50 for s Indi Ba and 880-940 K and 5.25 for s Indi 
Bb, both with a metallicity of [M/H]=-0.2. 

COND03 evolutionary model predictions for the masses are 
significantly inconsistent with the measured system mass if the 
young age range of 0.8-2.0 Gyr suggested by Lachaume et al. 
( 1999) is used. We find that a system age of 3.7^1.3 Gyr is neces- 
sary for the COND03 evolutionary models to be consistent with 
the measured system mass at the observed luminosities, and a re- 
view of the literature finds evidence supporting an age of ~5 Gyr 
for s Indi A. In the age range 3.7^1.3 Gyr, the COND03 mod- 
els predict effective temperatures in the range 1352-1385 K and 
976-101 1 K, for Ba and Bb, respectively. 

It is clear that there are several areas in which the atmo- 
spheric models currently do not reproduce observations and a 
more detailed analysis of these issues will be the subject of fu- 
ture work. They include the strength and shape of the wide ab- 
sorption by KI and Nal in the optical, the possible formation 
of feldspars in mid-late T dwarfs, and the reaction rates of CO 
and CH4. In addition, the spectral shape in the L-band caused 
by CH4 absorption is poorly reproduced, as is also the case for 
CH4 absorption at ~1.6/mi. The M-band spectra, although low 
resolution, also show that the atmospheric models significantly 
over-estimate the flux in this region. While the flux levels of the 
near-IR peaks can be reasonably reproduced, the level of absorp- 
tion between the peaks tends to be problematic. In particular, we 
find a feature at 1.35-1. 40 yum in both our object spectra which 
is not predicted in the atmospheric models. 

Neither source has detectable atomic Li I absorption at 
6708 A. The absence of lithium in the more massive compo- 
nent is consistent with the revised, higher age estimates coupled 
with its probable dynamical mass, while the lack of absorption in 
the cooler source is expected from its low effective temperature, 
where lithium is incorporated into molecules. 
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Although there is significant room for improvement in the 
atmospheric models, the current match to s Indi Ba and Bb is 
nevertheless impressive. When new data on methane opacities 
become available, we will be able to better reproduce the ob- 
served spectra and more reliably compare these spectral models 
to spectra of objects with less well-constrained physical param- 
eters. Additionally, when these updated atmospheric models are 
incorporated into the evolutionary models, a fully self-consistent 
comparison will be possible. 

Finally, when the individual dynamical masses become 
available and if we can obtain a reliable estimate of the age 
of this system, based on asteroseismological observations of the 
parent star s Indi A, then s Indi Ba and Bb will become invalu- 
able benchmark objects with a full set of physical parameters 
which newer models will have to reproduce, making them more 
reliable for analysing the properties of isolated ultra-cool field 
dwarfs. 

The predictions of the evolutionary models using luminosity 
and mass constraints are somewhat different to the derived ef- 
fective temperature and surface gravity from fitting atmospheric 
models to observed spectra. These differences may be resolved 
when the newer atmosphere models are incorporated into the 
evolutionary models. However, it seems that derivations of the 
mass of cool brown dwarfs are uncertain even where estimates 
of the effective temperature, surface gravity, and luminosity ex- 
ist. We therefore caution against the over-analysis of predicted 
brown dwarf masses at this time. 
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Appendix A: Image Fitting with Analytic Functions 

The 2-dimensional profile, or the point spread function (PSF), 
of ground-based optical/IR images is a superposition of several 
effects ( |Racine|1996| l, most importantly the spreading effect of 
the Earth's atmosphere caused by the mixing of air of different 
temperatures, leading to different refractive indices. In addition, 
there is the contribution of the optics and detectors used and any 



unintentional telescope motions, for example, due to imperfect 
tracking, which can cause elliptical profiles. 

Various authors have attempted to produce an analytical 
function which represents the shape of the PSF (eg. |King|197"T) 
Bendinelli et al. 1987 1. In a study of the two-dimensional profile 
of stellar images on photographic plates, Moffat] ( |1969| l found 
that a Gaussian profile, usually assumed to be a good match to 
the seeing, under-estimated the flux from the star at large ra- 
dial distances. He proposed an analytical profile (now termed the 
Moffat profile) to be a better match to photographic stellar im- 
ages. Similarly, King| ( 1971] ) found that a Gaussian profile was 
not sufficient to match observations, proposing instead a profile 
composed of a Gaussian core, falling to an exponential which 
tails to a inverse-square aureole. Franz] ( 1973[ ) presented a fur- 
ther analytical representation of the PSF which is often referred 
to as a modified Lorentzian. This was later successfully applied 
to CCD images of stellar profiles by Diego] (jT985j. 

We used the Levenberg-Marquardt technique to extract the 
best-fit parameters for each of three different analytical PSF 
models: Gaussian, elliptical Moffat, and elliptical modified 
Lorentzian profiles, with the aim of extracting the ratio of the 
fluxes of the two objects, thus allowing individual magnitudes to 
be found. The pixelation of the profile was accounted for in the 
fitting routine by pixelating the model values to the same resolu- 
tion as the data. The Moffat profile used is an elliptical version 
of the profile proposed by Moffat ( |1969| l: 
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where A is the peak amplitude, (xo,yo) the central co- 
ordinates, and a\, ai, and 6 define the semi -major and semi- 
minor axes, and the position angle of the ellipse. The elliptical 
modified Lorentzian used followed that of Diego ( 1985 1: 



L(x,y) 



where a x and ct2 are scale factors which relate to the full- 
width at half-maximum (Diego's RX and RY), and and a\ 
are similar scale factors which allow the exponent of By to vary 
with position (Diego's PRX and PRY). 

In Table [T] we showed the flux ratio and uncertainty derived 
in each band for both objects along with the central wavelengths 
and widths of the observed filters. The uncertainty on the fitted 
amplitudes for each image was found from the covariance matrix 
calculated by the Levenberg-Marquardt fit, which is a reasonable 
estimate of the standard error when, as in this case, the x 2 func- 
tion is quadratic. However, the uncertainties on the reported flux 
ratios are the standard errors of the fits for all images for each 
filter, which were comparable to the uncertainties from the co- 
variance matrix. The uncertainties on the flux ratios are not the 
dominant uncertainties in the final photometry of the two objects 
and so the anti-correlation is negligible. 
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Fig. A.l. An ISAAC /-band image of the s Indi Ba, Bb system 
(top left) and the residuals of the Gaussian (top right), modified 
Lorentzian (bottom left), and Moffat profile fits (bottom right). 
Each image is a 20 x 20 pixel (2.96" x 2.96") sub-section of the 
full image. North is up, East left. The scale applies to the residual 
images and shows the flux as a fraction of the observed peak 
pixel flux. It is apparent that the Gaussian profile is the worst fit, 
while the Moffat profile is somewhat better than the Lorentzian. 



The residuals of the best-fit profiles for each of the three 
functions are shown in Fig. A.l for an example observation and 
demonstrate the suitability of each profile to match the data. The 
Gaussian is seen, as expected, to provide the worst fit, under- 
estimating the peak and wings while over-estimating in between 
as it attempts to match the entire profile. If this profile were used 
to fit the relative fluxes of isolated stars in an image, one would 
expect it to return reasonable results, but the inability of the 
Gaussian to determine the flux in the wings of an object which 
blends with the profile of another makes it ineffective at extract- 
ing relative photometry of blended objects. The residuals of the 
modified Lorentzian profile however, show much reduced levels, 
with the Moffat profile residuals even lower. 

As confirmation of the validity of our PSF-fitting routine, 
it was employed to fit our optical images which we had previ- 
ously fit with the DAOPHOT/IRAF algorithm of |Stetson| ( [T987 1. 
DAOPHOT uses an analytical profile to model the core of the 
PSF and builds up an empirical determination of the wings from 
the chosen model stars. Although not completely independent 
of our method, in that it too uses an analytical profile, it is the 
matching of the model PSF to the wings that allows DAOPHOT 
to accurately extract the flux of blended objects. We fit 5 images 
in the /- and z-band and 12 images in the /?-band. Our fitting 
routine finds flux ratios of 4.71+0.08, 4.99+0.05, and 3.78+0.08 
for the R-, I-, and z-bands respectively, while with DAOPHOT 
we find flux ratios of 4.81+0.07, 5.04+0.05, and 3.85+0.03. The 
results of both PSF-fitting routines agree within the uncertain- 
ties in all cases, confirming the validity of our wholly analytical 
routine. Although the fitted flux ratio from DAOPHOT is larger 




Fig. B.l. The top panel shows the spectrum of s Indi Ba (upper 
line) and Bb in the region 0.975-1 .022 /urn where absorption by 
FeH and CrH is present. The second panel from the top shows 
an image of the observed spectrum (e Indi Ba is the lower spec- 
trum) sloping across the detector. The third panel from the top 
shows the corresponding fitted double Gaussian profile and the 
bottom panel shows the fitted double Moffat profile. The vertical 
striping seen in the fits is an effect of utilising the signal of all 
the pixels in the spatial direction to increase the signal-to-noise 
and so detect fine spectral features. 



than that from our fitting routine in these three examples, this is 
not true for all fitted images. 



Appendix B: Spectral Fitting Routine 

As with the imaging, the spectra of the two brown dwarfs were 
partially blended in the spatial direction and so a bespoke fitting 
algorithm was implemented to extract the individual spectra. The 
process was similar to the image fitting routine. We iteratively fit 
the parameters of a double Gaussian profile to the spatial direc- 
tion at each column along the spectra. Initially all parameters are 
free, but after the first fit we calculate the mean of the separations 
of the two peaks at each wavelength step, weighted by the total 
flux at each wavelength, from the ~1000 profiles in each image 
and refit the profiles with the separation fixed. We then do the 
same to derive a global fit to the spatial FWHM which we take 
to be constant with wavelength over each spectral window. The 
absolute position of the peaks at each wavelength was then con- 
strained by a trace of the spectrum across the detector so that 
an accurate centre would be found even for wavelengths with 
little signal-to-noise. With these constraints, the remaining pa- 
rameters were refit and the amplitudes were used along with the 
wavelength calibration to construct the individual wavelength- 
calibrated spectra. 

We found that a double Gaussian profile provided the best 
fit to our spectroscopic data, unlike in our broadband imaging 
where a Moffat profile was preferred. Figure B. 1 shows an image 
of the spectrum of s Indi Ba, Bb in the region 0.975-1. 022 /mi, 
along with fitted spectra, one using Gaussian profiles, the other 
Moffat profiles. The apparent vertical striping in the fitted spec- 
tra are real spectral features revealed by using the entire profile 
to increase the signal-to-noise. 

The resulting residuals for each profile (Fig. |B.2[ ) show 
clearly that the Gaussian profile is the better match. The Moffat 
profile under-estimates the peak of the profile while over- 
estimating the wings, consequently the contribution of the flux 
from one object to the other is not well determined. While an 
analytical Gaussian profile should not in principle be an exact 
match to the observed stellar spectral profile, the residuals are 



Robert R. King et al.: s Indi Ba, Bb: a detailed study 



27 




Fig. B.2. Residuals of the Gaussian (top) and Moffat profile fits 



to the spectral region shown in FigB.l The ill-fitting wings of 



the Moffat profile suggest that it is not a good approximation 
for the spatial profile in our spectroscopic images, while the 
Gaussian profile residuals are almost consistent with the back- 
ground noise. 
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Fig. C.l. The JHK band-passes of the MKO filter-set (red lines) 
and the ISAAC JHK$ filter-set (blue lines). The top-most line 
(thin, grey line) traces a typical atmosphere at Paranal showing 
the deep absorption bands between the filters. Also shown is the 
spectrum of a G8V star (thick grey line) from the Pickles spec- 
tral library (Pickles 1998) along with a smoothed spectrum of 
b Indi Ba (thick, black line). This highlights the difference be- 
tween the fraction of each object's flux in the regions of high at- 
mospheric absorption. Importantly, we see that the MKO filters 
mostly avoid these regions, while the ISAAC filters, especially 
the /-band, extend into them. Photometry using filters which ex- 
tend into these regions are susceptible to variations in the atmo- 
sphere and mismatches between target and standard star spectra. 



almost consistent with the background noise. It may have been 
expected that the same profile would be preferred for the spec- 
troscopy and broadband imaging, but if we consider imaging 
to be the integral of many Gaussian profiles with the FWHM 
varying as a function of wavelength, then the result would not 
be truly Gaussian. Nevertheless, the residual images and spec- 
tra clearly favour different profiles. The resulting full resolution 
spectra are presented in Figs.B]47] 



Appendix C: Photometric Calibration 

To allow a meaningful comparison of our sources with other ob- 
served brown dwarfs, it is necessary to place the magnitudes 
on a common photometric system due to the differences be- 
tween filter systems. In the near- to thermal-IR, the Mauna Kea 
Consortium filter-set ( |Tokunaga et al.||2002|l has been c hosen 
for this purpose by several groups (Stephens & Leggett 2004; 
|Golimowski et al.|2004||Hewett et al.|2006| l since these filters do 



not extend into the water absorption bands of the atmosphere and 
so avoid the large differences in relative spectral response be- 
tween sites and from varying atmospheric conditions. The differ- 
ences between the MKO JHK filters and the ISAAC JHK S fil- 



ters used for our observations are shown in Fig. C. 1 Many obser- 



vations have also been reported in the 2MASS system (Carpenter 



2001 ) due to the large number of ultra-cool dwarfs found in the 
2MASS database, although this system still encroaches on the 
regions of telluric absorption. We therefore present our near-IR 
photometry of s Indi Ba, Bb in both MKO (Table[2} and 2MASS 
systems (Table |4jl to allow easy comparison with other data. In 
the optical, we observed e Indi Ba, Bb in the FORS2 Bessell V 
and /, the FORS2 R special, and the FORS2 Gunn z filters and 
report our photometry in the FORS2 system. 

As none of our standard stars were of similar colour to our 
targets, when transforming our photometry into standard sys- 
tems, we could not use standard colour equations, nor could we 
ignore colour terms due to the filter differences. This is clear 
from Fig. C.l where we see that a T dwarf has relatively little 
flux in the regions of high absorption in the Earth's atmosphere 
(due to the dominance of H2O in both cases), but the standard 
star (in this case a solar-type star) will lose a large fraction of its 
flux in these regions. As a result, slight differences in the near-IR 
profiles of different systems can result in magnitude differences 
as large as 0.1 m (cf. |Stephens & Leggett|2"004) l. 

In the near-IR, photometric calibration used the solar type 
star, S234-E ( |Persson et al.||T998l >, as our standard star in the 
JHKs bands with the magnitudes known in the LCO (Las 
Campanas Observatory) system. HD205772 (A3) and HR8042 
(G3IV) were used in the L-band and Mnb -bands, respectively. 
The L-band standard star magnitude was known in the old 
UKIRT L-band filter and was assumed to be unchanged on trans- 
formation to the MKO L'-band as for other A-type stars (see 
Fig. 2 of Leggett et al. 2003 and the UKIRT photometric cali- 
bration web-page^]), with the spread in the L/L' magnitude dif- 
ferences of A stars being used as an estimate of the uncertainty 
on the transformation. For the Mnb -band, the standard star mag- 
nitude was known in the ESO system of van der Bliek et al. 
( fT9%l l. 

With the flux calibrated spectra of s Indi Ba and Bb, we 
extracted synthetic photometry in the MKO JHKL'M' filters 
and the 2MASS JHK$ filters by convolving the flux-calibrated 
spectra with the appropriate filter profiles and atmosphere for 
each site. In deriving our synthetic magnitudes, we used the 
ISAAC filters convolved with the model atmosphere given in the 
ISAAC user manuaQfor typical conditions over Paranal. For the 
2MASS magnitudes, we used the relative spectral responses of 
Cohe n et aL] ( [20 03 ) and applied the zero-point offsets of +0.001, 
-0.019, +0.017 for the J, H, and K s -bands respectively. The 
MKO filters used were convolved with the 1 .2 mm PWV (precip- 
itable water vapour) ATRAN model atmosphere of Lord ( 1992 ), 
as used by Stephens & Leggett (2004). We did not include the 
quantum efficiency of the detector nor the transmission profiles 
of any other optical elements. These were assumed to be prac- 
tically flat across each near-IR filter as discussed in Stephens & 
Leggett (2004). However, although the reflectivity of aluminium 
(used to coat the VLT mirrors) is relatively constant across the 
near-IR^] we note that it is strongly wavelength-dependent in the 



4 http://www.jach.hawaii.edu/UKIRT/astronomy/calib/phot_cal- 
/ukirt_stds.html 

5 http://www.eso.org/sci/facilities/paranal/instruments/isaac/doc/ 

6 http://www.gemini.edu/files/docman/press_releases/pr2004- 
5/images/comparison_AgAl_02.GIF 
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optical regime. We therefore did not attempt to derive our optical 
photometry in the same manner. 

That said, the precise optical filters used significantly affect 
the derived photometry as the flux of a T dwarf rises by more 
than three orders of magnitude over the range 0.6-1 .0 jum. Since 
our observed spectra do not cover the full V- and R -bands, we 
have not attempted to transform our optical photometry via syn- 
thetic photometry. Instead, we assumed that the magnitudes of 
our st andard stars, a s given by Landolt (1992} (VRI in the sys- 
tem of |Bessell| ( | 1 990] l) and SDSS (z), were the same as that in the 
FORS2 system as would be the case for an AO star. The spread in 
derived zero-points for the different standard stars were smaller 
than the other uncertainties which leads us to believe that the 
precise spectral type of the standard stars did not significantly 
affect the derived photometry. However, our optical photometry 
of s Indi Ba, Bb is not what would be measured through the stan- 
dard Bessell filters and so is not readily compared with other T 
dwarf observations. 



